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Confirming previous heuristic analyses a la Belinskii-Khalatnikov-Lifshitz, it is 
rigorously proven that certain "subcritical" Einstein-matter systems exhibit a mono- 
tone, generalized Kasner behaviour in the vicinity of a spacelike singularity. The 
D— dimensional coupled Einstein-dilaton-p-form system is subcritical if the dilaton 
couplings of the p-forms belong to some dimension dependent open neighbourhood 
of zero [|[], while pure gravity is subcritical if D > 11 O]. Our proof relies, like 
the recent theorem |T^] dealing with the (always subcritical [jl4j) Einstein-dilaton 
system, on the use of Fuchsian techniques, which enable one to construct local, an- 
alytic solutions to the full set of equations of motion. The solutions constructed are 
"general" in the sense that they depend on the maximal expected number of free 
functions. 



1 Introduction 



1.1 The problem 

In recent papers [p], @, 0], the dynamics of the coupled Einstein-dilaton-p-form system in 
D spacetime dimensions, with action (in units where 8irG = 1), 

% a/3 ,0,4?... 7n .] = 5 B [^]+^[^,0] + ^^[^,0,^ 1 ).. 7n .]+ "more", (1.1) 
S^] = ^jR^9d D x, (1.2) 
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^[^,0] = --Jd^d^V^gd n x, (1.3) 

^[^,0,^... 7nj ] = -_^|fW.^ (1.4) 

was investigated in the vicinity of a spacelike ( "cosmological" ) singularity along the lines 
initiated by Belinskii, Khalatnikov and Lifshitz (BKL) ||. In ( |1.1[ ), g a p is the spacetime 
metric, 6 is a massless scalar field known as the "dilaton", while the A%) „, are a col- 

lection of A; exterior form gauge fields (j = 1, ■ ■ ■ , k), with exponential couplings to the 
dilaton, each coupling being characterized by an individual constant Xj ( "dilaton coupling 
constant"). The F^'s are the exterior derivatives F^> = dA^\ whereas "more" stands for 
possible coupling terms among the p-forms which can be either of the Yang-Mills type (1- 
forms), Chern-Simons type || or Chapline-Manton type || [TJ. The degrees of the p-forms 
are restricted to be smaller than or equal to D — 2 since a (D — l)-form (or D-form) gauge 
field carries no local degree of freedom. In particular, scalars (nj = 0) are allowed among 
the A^^s but we then require that the corresponding dilaton coupling Xj be non-zero, so 
that there is only one "dilaton". Similarly, we require Xj ^ for the (D — 2)-forms (if any), 
since these are "dual"[| to scalars. This restriction to a single dilaton is mostly done for 
notational convenience: if there were other dilatons among the 0-forms, then, these must 
be explicitly treated on the same footing as and separated off from the p-forms because 
they play a distinct role. In particular, they would appear explicitly in the generalized 
Kasner conditions given below and in the determination of what we call the subcritical 
domain. The discussion would proceed otherwise in the same qualitative way. 

The main motivation for studying actions of the class ( |1.1| ) is that these arise as bosonic 
sectors of supergravity theories related to superstring or M-theory. In fact, in view of 
various no-go theorems, p-form gauge fields appear to be the only massless, higher spin 
fields that can be consistently coupled to gravity. Furthermore, there can be only one type 
of graviton || . With this observation in mind, the action ( [Lip is actually quite general. 
The only restriction concerns the scalar sector: we assume the coupling to the dilaton 
to be exponential because this corresponds to the tree-level couplings of the dilaton field 
of string theory. Note, however, that string-loop effects are expected to generate more 

1 We recall that the Hodge duality between a (nj + l)-form and a (D — rij — l)-form allows one to 
replace (locally) a n^-form potential by a (D — nj — 2)-form potential ) (with dilaton coupling 

K -V- 
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general couplings exp(A0) — > B(<p) which can exhibit interesting "attractor" behaviours 
||. We also restrict ourselves by not including scalar potentials; see, however, the end of 
the article for some remarks on the addition of a potential for the dilaton, which can be 
treated by our methods. 

Two possible general, "competing" behaviours of the fields in the vicinity of the space- 
like singularity have been identified^: 

1. The simplest is the "generalized Kasner behaviour" , in which the spatial scale factors 
and the field exp(0) behave at each spatial point in a monotone, power-law fashion 
in terms of the proper time as one approaches the singularity, while the effect of the 
p-form fields A^'s on the evolution of g^ y and <p can be asymptotically neglected. 
In that regime the spatial curvature terms can be also neglected with respect to the 
leading order part of the extrinsic curvature terms. In other words, as emphasized 
by BKL, time derivatives asymptotically dominate over space derivatives so that 
one sometimes uses the terminology "velocity-dominated" behaviour [PHI , instead of 
"generalized Kasner behaviour". We shall use both terminologies indifferently in 
this paper, recalling that in the presence of p-forms, which act as potentials for the 
evolution of the spatial metric and the dilaton (as do the spatial curvature terms), 
"velocity-dominance" means not only that the spatial curvature terms can be ne- 
glected, but also that the p-forms can be neglected in the Einstein-dilaton evolution 
equations.^ 



2. The second regime, known as "oscillatory" [||, or "generalized mixmaster" \Y1\ be 



haviour, is more complicated. It can be described as the succession of an infinite 
number of increasingly shorter Kasner regimes as one goes to the singularity, one 
following the other according to a well-defined "collision" law. This asymptotic evo- 
lution is presumably strongly chaotic. It is expected that, at each spatial point, the 
scale factors of a general inhomogeneous solution essentially behave as in certain ho- 
mogeneous models. For instance, for D = 4 pure gravity this guiding homogeneous 
model is the Bianchi IX model [EJ, while for D — 11 supergravity it is its naive 



one-dimensional reduction involving space-independent metric and three- forms 0. 

Whether it is the first or the second behaviour that is relevant depends on: (i) the 
spacetime dimension D, (ii) the field content (presence or absence of the dilaton, types of 
p-forms), and (iii) the values of the various dilaton couplings Xj. Previous work reached 
the following conclusions: 

• The oscillatory behaviour is general for pure gravity in spacetime dimension 4 
in fact, in all spacetime dimensions 4 < D < 10, but is replaced by a Kasner-like 
behaviour in spacetime dimensions D > 11 [|13|]. (The sense in which we use "general" 
will be made precise below.) 



2 For a recent extension of these ideas to the brane- worlds scenarios, see p"o[ . 

3 The Kasner solution is generalized in two ways: first, the original Kasner exponents include a dilaton 
exponent (if there is a dilaton), which appears in the Kasner conditions; second, the exponents are not 
assumed to be constant in space. We shall shorten "exhibits generalized Kasner behaviour" to Kasner-like. 
We stress that we do not use this term to indicate that the solution becomes asymptotically homogeneous 
in space. 
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The Kasner-like behaviour is general for the gravity-dilaton system in all spacetime 



dimensions D > 3 (see [|1J, fij| for D = 4). 



The oscillatory behaviour is general for gravity coupled to p-forms, in absence of a 
dilaton or of a dual (D — 2)-form (0 < p < D — 2) 0. In particular, the bosonic 
sector of 11-dimensional supergravity is oscillatory ]]J. Particular instances of this 
case have been studied in [|6|, |17|, [L8[ . 



• The case of the gravity-dilaton-p-form system is more complicated to discuss because 
its behaviour depends on a combination of several factors, namely the dimension 
D, the menu of p-forms, and the numerical values of the dilaton couplings. For 
a given D and a given menu of p-forms there exists a "subcritical" domain T> (an 
open neighbourhood of the origin Xj = for all j's) such that: (i) when the Xj 
belong to D the general behaviour is Kasner-like, but (ii) when the Xj do not belong 
to T> the behaviour is oscillatory. Note that T> is open. Indeed, the behaviour is 
oscillatory when the Xj are on the boundary of V, as happens for instance for the 
low-energy bosonic sectors of type I or heterotic superstrings ||]. For a single p-form, 
the subcritical domain V takes the simple form \Xj\ < X C j, where A^ depends on the 
form-degree and the spacetime dimension. (A^ can be infinite.) For a collection of 
p-forms, T> is more complicated and not just given by the Cartesian product of the 
subcritical intervals associated with each individual p-form. 

The above statements were derived by adopting a line of thought analogous to that 
followed by BKL. Now, as understood by BKL themselves, these arguments, although 
quite convincing, are somewhat heuristic. It is true that the original arguments have 



received since then a considerable amount of both numerical and analytical support |19 



20] , |21| , |22| , |23| . Yet, they still await a complete proof. One notable exception is the four- 



dimensional gravity-dilaton system, which has been rigorously demonstrated in |L5| to be 
indeed Kasner-like, confirming the original analysis []nj. Using Fuchsian techniques, the 



authors of [|15[| have proven the existence of a local (analytic) Kasner-like solution to the 
Einstein-dilaton equations in four dimensions that contains as many arbitrary, physically 
relevant functions of space as there are local degrees of freedom, namely 6 (counting q 
and q independently). To our knowledge, this was the first construction, in a rigorous 
mathematical sense, of a general singular solution for a coupled Einstein-matter system. 
Note in this respect several previous works in which formal solutions had been constructed 
near (Kasner-like) cosmological singularities by explicit perturbative methods, to all orders 
of perturbation theory p4j, ^5f . 

The situation concerning the more complicated (and in some sense more interesting) 
generalized mixmaster regime is unfortunately - and perhaps not surprisingly - not so well 
developed. Rigorous results are scarce (note p6[] ) and even in the case of the spatially 



homogeneous Bianchi IX model only partial results exist in the literature [27 



The purpose of this paper is to extend the Fuchsian approach of [[HJ to the more 



complicated class of models described by the action (1J.) and to prove that those among 
the above models that were predicted in |TB|, [TJ, |2j to be Kasner-like are indeed so. This 
provides many new instances where one can rigorously construct a general singular solution 



3 



for a coupled Einstein-matter (or pure Einstein, in D > 11) system. In fact, our (Fuchsian- 
system-based) results prove that the formal perturbative solutions that can be explicitly 
built for these models do converge to exact solutions. This provides a further confirmation 
of the general validity of the BKL ideas. We shall also explicitly determine the subcritical 
domain T> for a few illustrative models. For all the relevant systems, we construct local 
(near the singularity) analytic solutions, which are "general" in the sense that they contain 
the right number of freely adjustable arbitrary functions of space (in particular, these 
solutions have generically no isometries), and which exhibit the generalized (monotone) 
Kasner time dependence. 



1.2 Strategy and outline of the paper 



Our approach is the same as in [pL5| , and results from that work will be used frequently 
here without restating the arguments. Here is an outline of the key steps. 

A d + 1 decomposition is used, for d spatial dimensions, d = D — 1. A Gaussian 
time coordinate, t, is chosen such that the singularity occurs at t — 0. The first step 
in the argument consists of identifying the leading terms for all the variables. This is 
accomplished by writing down a set of evolution equations which is obtained by truncating 
the full evolution equations, and then solving this simpler set of evolution equations. This 
simpler evolution system is called the Kasner-likeQ evolution system (or, alternatively, the 
velocity-dominated system). It is a system of ordinary differential equations with respect 
to time (one at each spatial point) which coincides with the system that arises when 
investigating metric-dilaton solutions that depend only on time. The precise truncation 
rules are given in subsection |2.2| below. The second step is to write down constraint 
equations for the Kasner-like system (called "velocity-dominated" constraints) and to show 
that these constraints propagate, i.e., that if they are satisfied by a solution to the Kasner- 
like evolution equations at some time t > 0, then they are satisfied for all time t > 0. In the 
set of Kasner-like solutions, one expects that there is a subset, denoted by V, of solutions 
which have the property of being asymptotic to solutions of the complete Einstein-dilaton- 
p-form equations as t — > 0, i.e as one goes to the singularity. This subset is characterized 
by inequalities on some of the initial data, which, however, are not always consistent. The 
existence of a non-empty V requires the dilaton couplings to belong to some range, the 
"subcritical range". When V is non-empty and open, the solutions in V involve as many 
arbitrary functions of space as a "general solution" of the full Einstein equations should. 
On the other hand, if V is empty the construction given in this paper breaks down and 
the dynamical system is expected to be not Kasner-like but rather oscillatory. 

To show that indeed, the solutions in V (when it is non-empty) are asymptotic to 
true solutions, the third step is to identify decaying quantities such that these decaying 
quantities along with the leading terms mentioned above uniquely determine the variables, 
and to write down a Fuchsian system for the decaying quantities which is equivalent to 
the Einstein-matter evolution system. As the use of Fuchsian systems is central to our 
work let us briefly recall what a Fuchsian system is and how such a system is related to 



4 Note that we use the terms "Kasner-like solutions" to label both exact solutions of the truncated 
system and solutions of the full system that are asymptotic to such solutions. Which meaning is relevant 
should be clear from the context. 
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familiar iterative methods. For a more detailed introduction to Fuchsian techniques see 
fl5| , |28| , p9| |30] and references therein. Note that we shall everywhere restrict ourselves to 
the analytic case. We expect that our results extend to the C°° case, but it is a non-trivial 
task to prove that they do. 

The general form of a Fuchsian system for a vector-valued unknown function u is 

t d t u + A(x) u = f(t, x, u, u x ), (1.5) 

where the matrix A(x) is required to satisfy some positivity condition (see below), while 
the "source term" / on the right hand side is required to be "regular." (See [|15] for 
precise criteria allowing one to check when the positivity assumption on A(x) is satisfied 
and when / is regular.) A key point is that / is required to be bounded by terms of 
order 0(t 5 ) (with t -> 0,5 > 0) as soon as u and their space derivatives u x are in a 
bounded set (a simple, concrete example of a source term satisfying this condition is 
f = t Sl + t S2 u + t Si u x , with <5j's larger than 5). A convenient form of positivity condition 
to be satisfied by the matrix A(x) is that the operator norm of be bounded when 
< t < 1 (and when x varies in any open set). Essentially this condition restricts the 
eigenvalues of the matrix A(x) to have positive real parts. The basic property of Fuchsian 
systems that we shall use is that there is a unique solution to the Fuchsian equation which 
vanishes as t tends to zero [28|. One can understand this theorem as a mathematically 



rigorous version of the recursive method for solving the equation ( |1.5|) . Indeed, when 



confronted with equation ( |1.5| ), it is natural to construct a solution by an iterative process, 
starting with the zeroth order approximation uq = (which is the unique solution of 
( |1.5[ ) with / = that tends to zero as t — > 0), and solving a sequence of equations of 
the form td t u^ n ' + A{x)u^ n ' = f(t,x,u^ n ~ 1 \u x n ^). At each step in this iterative process 
the source term is a known function which essentially behaves (modulo logarithms) like a 
sum of powers of t (with space-dependent coefficients). The crucial step in the iteration 
is then to solve equations of the type tdtu + A(x)u = C(x)t &<yX \ The positivity condition 
on A(x) guarantees the absence of homogeneous solutions remaining bounded as t — > 0, 
and ensures the absence of "small denominators" in the (unique bounded) inhomogeneous 
solution generated by each partial source term: Ui n h om = (5 + A)~ 1 Ct s . (See, e.g., [E| for 
a concrete iterative construction of a Kasner-like solution and the proof that it extends to 
all orders.) This link between Fuchsian systems and "good systems" that can be solved to 
all orders in a formal iteration makes it a priori probable that all cases which the heuristic 
approach a la BKL has shown to be asymptotic to a Kasner-like solution (by checking 
that the leading "post-Kasner" contribution is asymptotically sub-dominant) can be cast 
in a Fuchsian form. The main technical burden of the present work will indeed be to 
show in detail how this can be carried out for the evolution systems corresponding to all 
the sub-critical {i.e., non-oscillatory) Einstein-matter systems. Our Fuchsian formulation 
proves that (in the analytic case) the formal all-orders iterative solutions for the models we 
consider do actually converge to the unique, exact solution having a given leading Kasner 
asymptotic behaviour as t — ► 0. 

Finally, the fourth step of our strategy is to prove that the constructed solution does 
satisfy also all the Einstein and Gauss-like constraints so that it is a solution of the full 
set of Einstein-matter equations. We shall deal successively with the matter (Gauss-like) 
constraints, and the Einstein constraints. 
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Our paper is organized as follows. In section ^], we first consider the paradigmatic 
example of gravity coupled to a massless scalar field and to a Maxwell field in 4 spacetime 



dimensions. The action (1.1) reads in this case 



% a/3 ,0,Ay] = lJ{R-d^d^-^F,„F^e x *}V=gd 4 x. (1.6) 

For this simple example, we shall explicitly determine the subcritical domain V, i.e., the 
critical value A c such that the system is Kasner-like when — A c < A < A c . Because this case 
is exemplary of the general situation, while still being technically rather simple to handle, 
we shall describe in some detail the explicit steps of the Fuchsian approach. 

In section [J[ vacuum solutions governed by the pure Einstein action ( |1.2| ) with D > 1 1 
are considered. This system was argued in [|13| to be Kasner-like and we show here how this 
rigorously follows from the Fuchsian approach. Note that, contrary to what happens when 
a dilaton is present, Fuchsian techniques apply here even though not all Kasner exponents 
can be positive. 

In sections [| - 18|, the results of the previous sections are generalized to the wider class 
of systems (|1.1| ). First, in section [|, to solutions of Einstein's equation with spacetime 
dimension D > 3 and a matter source consisting of a massless scalar field, governed by the 
action SE[ga/3]+S^[g a /3, 4>}- This is the generalization to any D > 3 of the case D = 4 treated 
in |jT5[ . In section we turn to the general situation described by the action (|1.1|), without, 
however, including the additional terms represented there by "more" . We then give some 
general rules for computing the subcritical domain of the dilaton couplings guaranteeing 
velocity-dominance (section ^). The inclusion of interaction terms is considered in the last 
sections. It is shown that they do not affect the asymptotic analysis. This is done first for 
the Chern-Simons and Chapline-Manton interactions in section [!], and next, in section |8], 
for the Yang-Mills couplings (for some gauge group G), for which the action reads 

S[g aP , 0, A 1 ] = ~J{R- d,<P - -i^, • F^e x ^}^g d D x. (1.7) 

Here the dot product, F ■ F, is a time-independent, Ad-invariant, non-degenerate scalar 
product on the Lie algebra of G (such a scalar product exists if the algebra is compact, or 
semi-simple). Contrary to what is done in sections |5| and |7|, we must work now with the 
vector potential (and not just with the field strength), since it appears explicitly in the 
coupling terms. 

In section [| we show that self-interactions of a rather general type for the scalar field 
can be included without changing the asymptotics of the solutions. Explicitly, we add a 
(nonlinear) potential term, 

Sm\9ap,<l>\ = -J V(4>)^gd D x, (1.8) 

to the action ( |1 . 1| ) , where V(<p) must fulfill some assumptions given in section ^. V(<p) 
may, for example, be an exponential function of </>, a constant, or a suitable power of <fi. 
Similar forms for V(<p) were considered with D = 4 in EI]. 



Finally, in section [H| we state two theorems that summarize the main results of the 
paper and give concluding remarks. 
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1.3 On the generality of our construction 

As we shall see the number of arbitrary functions contained in solutions to the velocity- 
dominated constraint equations is equal to the number of arbitrary functions for solutions 
to the Einstein-matter constraints. In this function-counting sense, our construction de- 
scribes what is customarily called a "general" solution of the system. Intuitively speaking, 
our construction concerns some "open set" of the set of all solutions (indeed, the Kasner- 
like behaviour of the solution is unchanged under arbitrary, small perturbations of the 
initial data, because this simply amounts to changing the integration functions). Note 
that, in the physics literature, such a "general" solution is often referred to as being a 
"generic" solution. However, in the mathematics literature the word "generic" is restricted 
to describing either an open dense subset of the set of all solutions, or (when this can be 
defined) a subset of measure unity of the set of all solutions. In this work we shall stick to 
the mathematical terminology. We shall have nothing rigorous to say about whether our 
general solution is also generic. However, we wish to emphasize the following points. 

First, let us mention that the set V of solutions to the velocity-dominated equations 
that are asymptotic to solutions of the complete equations is not identical to the set U 
of all solutions to the velocity-dominated constraint equations. The subset V C U is de- 
fined by imposing some inequalities on the free data. These inequalities do not change 
the number of free functions. Therefore the solutions in V are still "general". One can 
wonder whether there could be a co-existing general behaviour, corresponding to initial 
data that do not fulfill the inequalities. For instance, could such "bad" initial data lead to 
a generalized mixmaster regime? This is a difficult question and we shall only summarize 
here what is the existing evidence. There are heuristic arguments, supported by numerical 



study, [14], |32], |33|, [34J that suggest that if one starts with initial data that do not fulfill the 
inequalities, one ends up, after a finite transient period (with a finite number of "collisions" 
with potential walls), with a solution that is asymptotically velocity-dominated, for which 
the inequalities are fulfilled almost everywhere. In that sense, the inequalities would not 
represent a real restriction since there is a dynamical mechanism that drives the solution to 
the regime where they are satisfied. For the subcritical values of the dilaton couplings that 
make the inequalities defining V consistent, there is thus no evidence for an alternative 
oscillatory regime corresponding to a different (open) region in the space of initial dataQ. 
It has indeed been shown that the inequalities defining V are no restriction in a large 



spatially homogeneous class f27|. Such rigorous results are, however, lacking in the inho- 



mogeneous case. In fact, an interesting subtlety might take place in the inhomogeneous 
case. The heuristic arguments and numerical studies of |33|, [34]] suggest the possibility that 
the mechanism driving the system to V may be suppressed at exceptional spatial points 
in general spacetimes, with the result that the asymptotic data at the exceptional spatial 
points are not consistent with the inequalities we assume and lead to so-called "spikes". 
This picture has been given a firm basis in a scalar field model with symmetries but 
the status of the spikes in a general context remains unclear. 

Finally, since we only deal with spacelike singularities, the classes of solutions we con- 



The oscillatory regime may however be present for peculiar initial data, presumably forming a set 
of zero measure. For instance, gravity + dilaton is generically Kasner-like, but exhibits an oscillatory 
behaviour for initial data with <fi = (in D < 11 spacetime dimensions). 
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sider do not contain all solutions governed by the action fll.lj ). Other types of singularities 
(e.g. timelike or null ones) are known to exist. Whether these other types of singularities 
are general is, however, an open question. 



1.4 Billiard picture 

At each spatial point, the solution of the coupled Einstein- matter system can be pictured, in 
the vicinity of a spacelike singularity, as a billiard motion in a region of hyperbolic space |]36| , 
37 , H|, 38fl . Hyperbolic billiards are chaotic when they have finite volume and non chaotic 
otherwise. In this latter case, the "billiard ball" generically escapes freely to infinity after 
a finite number of collisions with the bounding walls. Subcritical Einstein-matter systems 
define infinite-volume billiards. The velocity-dominated solutions correspond precisely to 
the last (as t — > 0) free motion (after all collisions have taken place), in which the billiard 
ball moves to infinity in hyperbolic space. 



1.5 Conventions 

We adopt a "mostly plus" signature ( — h + + . . .)■ The spacetime dimension is D = d + 1. 
Greek indices range from to d, while Latin indices G {1, . . . , d}. The spatial Ricci tensor 
is labeled R and the spacetime Ricci tensor is labeled ^R. Our curvature conventions 
are such that the Ricci tensor of a sphere is positive definite. Einstein's equations read 
G a/ 3 = T a p, where G a p = R a p — Rg a ^j2 denotes the Einstein tensor and T a/3 denotes the 
matter stress-energy tensor, T a p = — (2/ v /^)5S' matter /<5g Q/3 , an d units such that 8nG = 1. 
The spatial metric compatible covariant derivative is labeled V a and the spacetime metric 
compatible covariant derivative is labeled ( D > V a . The velocity-dominated metric compati- 
ble covariant derivative is labeled °V a . According to the context, g denotes the (positive) 
determinant of g a b in d + 1-decomposed expressions, and the (negative) determinant of g^ v 
in spacetime expressions. Whenever t s or t~ 5 appears, 5 denotes a strictly positive number, 
arbitrarily small. We use Einstein's summation convention on repeated tensor indices of 
different variances. (When the need arises to suspend the summation conventions for some 
non-tensorial indices, we shall explicitly mention it.) In expressions where there is a sum 
that the indices do not indicate, all sums in the expression are indicated explicitly by a 
summation symbol. Indices on the velocity-dominated metric and the velocity-dominated 
extrinsic curvature are raised and lowered with the velocity-dominated metric. 



1.6 d + 1 decomposition 

Consider a solution to the Einstein's equations following from (|1 . 1| ) , consisting of a Lorentz 
metric and matter fields on a .D- dimensional manifold M which is diffeomorphic to (0, T) x 
S for a d- dimensional manifold, E, such that the metric induced on each t = constant 
hypersurface is Riemannian, for t G (0,T). Here D is an integer strictly greater than 
two. Furthermore, consider ad+l decomposition of the Einstein tensor, Gap, and the 
stress-energy tensor, T a/ 3, with a Gaussian time coordinate, t G (0,T), and a local frame 
{e a } on E. Note that the frame e a = e l a (x)di is time-independent. The spacetime metric 
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reads ds 2 = —dt 2 + g ab (t, x)e a e b , where e a = e°;(x)dx l (with e°;e l b = 5%) is the co-frame. Let 
P = T 00 , j a = -T 0a and S ab = T ab . Define 

C = 2G 00 -2Too (1.9) 
= -k a b k b a + (trk) 2 + R-2p. 

C = is the Hamiltonian constraint. Similarly, C a = is the momentum constraint, where 

C a = -G 0a + T 0a (1.10) 
= V b k b a -V a (trk)-j a . 

In Gaussian coordinates, the relation between the metric and the extrinsic curvature is 

d t g a b = -2k ab . (1.11) 

The evolution equation for the extrinsic curvature is obtained by setting E a b = 0, with 

E\ = t D )R\-T\+ —L_ T8 a b (1.12) 

d t k\ = R a b + (tr k)k\-M a h . (1.13) 

Here 

M a b = S a b -^-^((trS)-p)5 a b . 

2 Scalar and Maxwell fields in four dimensions 
2.1 Equations of motion 

As said above, let us start by considering in detail, as archetypal system, the system defined 
by the action (|1.6|) , i.e., the spacetime dimension is D = 4 and the matter fields are a 
massless scalar field exponentially coupled to a Maxwell field, with the magnitude of the 
dilaton coupling constant smaller in magnitude than some positive real number determined 
below, < |A| < A c . The stress-energy tensor of the matter fields is 

T, v = (4) (4) SJ v cj> - l -g^ u W V a (4) V> + [F ua F v a - ^F^F^e^. 
The matter fields satisfy the following equations. 

( 4 >Vj 4 )V> = ^F^e^, 

{4) V M (F^e A *) = 0, 
(4) V [Q ^ 7] = 0. 
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The 3+1 decomposition of the stress-energy tensor is best expressed in terms of the electric 
spatial vector density £ a = yfg F 0a e x ^ and the magnetic antisymmetric spatial tensor F ab . 

2 g 2 

j a = -d t( pe a {4>) + ^-S b F ab , 

M\ = g ac e b ( ( j ) )e c ( ( t>)--{g bc S a S c -l-5 a b g ch £ c £ h }e~^ 

9 2 

+{9 ac 9 hl F ch F bl - -S^^F^F^. (2.1) 

The matter constraint equations are 

e a (S a ) + f b ba S a = (2.2) 
e [a (F bc] ) + tf ab F c]h = 0. (2.3) 

Here f£ b are the (time-independent) structure functions of the frame, [e a , e b ] = /„ 6 e c . The 
matter evolution equations are 

d^-itrk)^ = g^VaV^+^ga^S'e-^-^g^FaAe^, (2.4) 

d t 8 a = e b (^g ac g bh F ch e x *) + (f ib g ac +^f^ (2.5) 

d t F ab = -2e [a (^-g b]c S c e-^) + f c ab ^-g ch £ h e~ x t (2.6) 

2.2 Velocity-dominated evolution equations and solution 

The Kasner-like, or velocity-dominated, evolution equations are obtained from the full 
evolution equations by: (i) dropping the spatial derivatives from the right hand sides of 
( |1.13j ), (|2.4|), ( |2.5| ) and ( |2.6| ) (note that /^ fe -terms count as derivatives and that we keep 
the time derivatives of the magnetic field in (|2.6j ) even though F ab = d a A b — d b A a ); and 
(ii) dropping the p-form terms in both the Einstein and dilaton evolution equations. This 
is a general rule and yields in this case 



9 t °g a b 


— 2 °k ab , 


(2.7) 


d t °k a b 


= (tr °k)°k a b , 


(2.8) 


d?°<J)-(tr k)d t <J) 


= o, 


(2.9) 


d t °£ a 


= o, 


(2.10) 


d t °F ab 


= 0. 


(2.11) 



(As we shall see below, interaction terms of Yang-Mills or other types - if any - should 
also be dropped.) 

It is easy to find the general analytic solution of the evolution system (|2.7p - Q2.ll 



since the equations are the same as for "Bianchi type I" homogeneous models (one such set 
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of equations per spatial point). Taking the trace of ( |2.8| ) shows that — l/tr°fc = t + C(x). 
By a suitable redefinition of the time variable one can set C(x) to zero. Then (|2.8| ) shows 
that —t°k a b = K a b is a constant matrix (which must satisfy trK = K a a = 1, and be such 
that g ac (to)K c b is symmetric in a and b), 



°k a b (t) = -r 1 K a b . (2.12) 

Injecting this information into ( |2.7| ) leads to a linear evolution system for °g ab - td t °g a b = 
2 °g ac K c b , which is solved by exponentiation, 



°9ab(t) = °gac{to) 



The other evolution equations are also easy to solve, 



(2.13) 



°<f>(t) = Alnt + B, (2.14) 
°S a (t) = °S a (2.15) 
°F ab (t) = °F ab . (2.16) 

In ( |2.13| ) (t/to) 2K denotes the exponentiation of a matrix. Quantities on the left hand 
side of (|2.12| ) - (|2.16|) may be functions of both time and space, while all the time de- 
pendence of the right hand side is made explicit. For instance, ( |2.16| ) is saying that the 
spacetime dependence of the general magnetic field °F ab (t, x) (solution of the velocity- 
dominated evolution system) is reduced to a simple space dependence °F ab (x) (where °F ab 
is an antisymmetric spatial tensor). Let p a denote the eigenvalues of K a b , ordered such 
that pi < P2 < P3- Since trK = 1, we have the constraint 

Pi+p 2 +p 3 = l- (2.17) 



In the works of BKL the matrix solution ( |2.13|) is simplified by using a special frame {e a } 
with respect to which the matrices °g a b(to) and K a b are diagonal. However, as emphasized 
in [ l5[ , this choice can not necessarily be made analytically on neighbourhoods where the 
number of distinct eigenvalues of K a b is not constant. To obtain an analytic solution, while 
still controlling the relation of the solution to the eigenvalues of K a b , a special construction 
was introduced in |15]| . This construction is based on some (possibly small) neighbourhood 



U of an arbitrary spatial point x G £ and uses a set of auxiliary exponents q a (x). These 
auxiliary exponents remain numerically close to the exact "Kasner exponents" p a {x), are 
analytic and enable one to define an analytic frame (see below). To construct the auxiliary 
exponents q a (x) one distinguishes three cases: 

Case I (near isotropic): If all three eigenvalues are equal at xq, choose a number 
e > so that for x G U , max ajb \p a ( x ) ~ Pb(x)\ < e /2. In this case define q a = 1/3 on U , 
a = 1,2,3. 

Case II (near double eigenvalue): If the number of distinct eigenvalues at xq is two, 
choose e > so that for x G Uo, max a ^ \p a — p b \ > e/2, and \p a > — pv\ < e/2 for some 
pair, a', b', a' ^ b', shrinking Uq if necessary. Denote by p± the distinguished exponent not 
equal to p a >, p b '. In this case define = p± and q a i = q b > = (1 — <7_l)/2 on Uq. 
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Case III (near diagonalizable): If all eigenvalues are distinct at xq, choose e > so 
that for x G U , min a b \p a (x) — p b (x)\ > e/2, shrinking U if necessary. In this case define 

q a = p a on U . 

The frame {e a }, called the adapted frame, is required to be such that the related (time- 
dependent) frame {e a (t) = t~ Qa e a } is orthonormal with respect to the velocity-dominated 
metric at some time t > 0, i.e., such that °g a b(to) — tl qa 5 a b- (Here and in the rest of the 
paper, the Einstein summation convention does not apply to indices on q a and p a . These 
indices should be ignored when determining sums. Furthermore, quantities with a tilde 
will refer to the frame {e a (t)}.) 

In addition, in Case II it is required that e± be an eigenvector of K corresponding to 
q± and that e a >, e b > span the eigenspace of K corresponding to the eigenvalues p a ', py. In 
case III it is required that the e a be eigenvectors of K corresponding to the eigenvalues 
q a {= Pa)- In all cases it is required that {e a } be analytic. The auxiliary exponents, q a , 
are analytic, satisfy the Kasner relation Y^Qa = 1, are ordered (q± < q<i < q%), and satisfy 
Qi > Pi, <?3 < P3 and max a \q a - p a \ < e/2. If q a ^ q b , then °g ab , °g ab , °g ab and °g ab all 
vanish, and the same is true with g replaced by k. 

Equations (|2.12|) - ( |2.16| ), with the form of g a b(t ) and K a b specialized as given just 



above, are the general analytic solution to the velocity-dominated evolution equations in 
the sense that any analytic solution to the velocity-dominated evolution equations takes 
this form near any xq G £ by choice of (global) time coordinate and (local) spatial frame. 

2.3 Velocity-dominated constraint equations 

When written in terms of the velocity-dominated variables, the velocity-dominated con- 
straints take the same form as the full constraint equations, except the Hamiltonian con- 
straint, which is obtained by dropping spatial gradients and electromagnetic contributions 
to the energy-density. This is a general rule, valid also for the more general models con- 
sidered below. Thus, if we define 



°j a = -^Ve a (V) + -^°^°^, 



g 

we get °C = and °C a = for the velocity-dominated constraints corresponding to the 
Hamiltonian and momentum constraints, with 

°C = -°k a b k b a + (ti°k) 2 -2°p, (2.18) 
°C a = °V b k b a -e a (tr k)- j a . (2.19) 

The velocity-dominated matter constraint equations read 

e a (°S a ) + f b ba °S a = 0, 
e[ a (°F bc] ) + f [ab °F c ] h = 0. 



12 



For the solution Q2.12| ) - ( |2.14j ) the velocity-dominated Hamiltonian constraint equation 
is equivalent to 

J2Pa 2 + A 2 = l. (2.20) 

The conditions ( ]2.17|) and ( |2.2(J| ) are the famous Kasner conditions when the dilaton is 
present. While p\ is necessarily non-positive when A = 0, this is no longer the case when 
the dilaton is nontrivial (A ^ 0): all p a 's can then be positive. This is the major feature 
associated with the presence of the dilaton, which turns the mixmaster behaviour of (4- 
dimensional) vacuum gravity into the velocity-dominated behaviour. We shall call (p a , A) 
the Kasner exponents (because they are the exponents of the proper time in the solution 
for the scale factors and exp0) and refer to (|2.17|) and fl2.20| ) as the Kasner conditions 



(note that A is often denoted p^ to emphasize its relation to the kinetic energy of 0, and 
its similarity with the other exponents). 
A straightforward calculation shows that 

<9 t °C-2(tr°A;) C = 0, (2.21) 
d t °C a -(tr°k)°C a = -~e a (°C). (2.22) 

Thus if the velocity-dominated Hamiltonian and momentum constraints are satisfied at 
some to > 0, then they are satisfied for all t > 0. Similarly, since °S a and F a j 1 are 
independent of time, if the matter constraints are satisfied at some time to > 0, then they 
are clearly satisfied for all time. 

2.4 Critical value of dilaton coupling A c 

Our ultimate goal is to show that the velocity-dominated solutions asymptotically approach 
(as t — > 0) solutions of the original system of equations. We shall prove that this is the 
case provided the Kasner exponents pi, A, subject to the Kasner conditions 

E^ 2 -(E^) 2 + ^ 2 = °> £p a = i ( 2 - 23 ) 

obey additional restrictions. These restrictions are inequalities on the Kasner exponents 
and read explicitly 

2pi - XA > 0, pi > 0, 2pi + XA > 0. (2.24) 

As explained in [Q], and rigorously checked below, these restrictions are necessary and suf- 
ficient to ensure that the terms that are dropped when replacing the full Einstein-dilaton- 
Maxwell equations by the velocity-dominated equations become indeed negligible as t —>■ 0. 
More precisely, the first condition (respectively the third) among ( |2.24| ) guarantees that 



one can neglect the electric (respectively, magnetic) part of the energy-momentum tensor 
of the electromagnetic field in the Einstein equations, whereas the condition pi > is 
necessary for the spatial curvature terms to be asymptotically negligible. The conditions 
( |2.24| ) define the set V of velocity-dominated solutions referred to in the introduction. 
It is clear that if |A| is small enough - in particular, if A = - the inequalities (|2.24|) 



can be fulfilled since the Kasner exponents can be all positive when the dilaton is included. 
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But if |A| is greater that some critical value A c , it is impossible to fulfill simultaneously 
the Kasner conditions (|2.23|) and the inequalities ( [2.241) , because one of the terms ±XA 



becomes more negative than 2pi is positive. In that case, the set V is empty and our 
construction breaks down. For |A| < A c , however, the set V is non-empty and, in fact, 
stable under small perturbations of the Kasner exponents since ( |2.24j ) defines an open 



region on the Kasner sphere. We determine in this subsection the critical value A c such 
that fl2.23|) and ( p.24| ) are compatible whenever |A| < A c . 



To that end, we follow the geometric approach of [|3], 3S]. In the 4-dimensional space of 
the Kasner exponents (p a , A), we consider the "wall chamber" W defined to be the conical 
domain where 

Pi<P2< Pa, 2p a -XA>0, Pa > 0, 2p a + \A > 0. (2.25) 
These inequalities are not all independent since the four conditions 

Pi<P2< Ps, 2pi - XA>0, 2ft + XA > (2.26) 
imply all others. The quadratic Kasner condition (|2.23| ) can be rewritten 



GVpV = 0, (jf) = {p a ,A) (2.27) 
where defines a metric in "Kasner-exponent space" 

dS 2 = G^dp»dp v = dpa 2 ~ (E d P«) 2 + ( dA T ( 2 - 28 ) 

The metric (|2.28|) has Minkowskian signature (— , +, +, +). An example of timelike direc- 
tion is given by p\ = p 2 = p 3 , A = 0. Inside or on the light cone, the function J^Pa does 
not vanish. The upper light cone (in the space of the Kasner exponents) is conventionally 
defined by (|2.27|) and the extra condition J^Pa > 0. It is clear from our discussion that the 



Kasner conditions ( |2.23| ) and the inequalities fl2.24Q are compatible if and only if there are 



lightlike directions in the interior of the wall chamber W (by rescaling p^ — > ap 11 , a > 0, 
one can always make ^ p a = 1 for any point in the interior of the wall chamber so that this 
condition does not bring a restriction). The problem amounts accordingly to determining 
the relative position of the light cone ( |2.27| ) and the wall chamber ( |2.26| ). 

This is most easily done by computing the edges of (|2.26|) , i.e., the one-dimensional 



intersections of three faces among the four faces (p.26|) of W. There are four of them: (i 



pi = p 2 = A = 0, p 3 = a; (ii) p x = A = 0; p 2 = p 3 = a; (iii) 2pi = 2p 2 = 2p 3 = XA = a; 
and (iv) 2pi = 2p 2 = 2p 3 = —XA = a, where in each case, a > is a parameter along the 
edge (a = being the origin). The vectors (A = 1, 2, 3, 4) along the edges corresponding 
to a = 1, namely (0, 0, 1, 0), (0, 1, 1, 0), (1/2, 1/2, 1/2, 1/A) and (1/2, 1/2, 1/2, -1/A) form 
a basis in Kasner-exponent space. Any vector t> M can thus be expanded along the e^, 
v ti = v A e ^ A _ \ point P in Kasner-exponent space is on or inside the wall chamber W if 
and only if its coordinates p A in this basis fulfill p A > with P inside when p A > for 
all A's. Thus, if all the edge vectors are timelike or lightlike, the Kasner conditions 
are incompatible with the inequalities fl2.24p since any linear combination of causal vectors 



with non-negative coefficients is on or inside the forward light cone (the e^'s are future- 
directed since p\ + p 2 + p 3 > for all of them). If, however, one (or more) of the edge 
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vectors lies outside the light cone, then, the Kasner conditions and the inequalities ( |2.24| ) 
are compatible. The nature of some of the edge vectors depends on the value of the dilaton 
coupling A: while the first one is always lightlike and the second one always timelike, the 
squared norm of the last two is -3/2 + 1/A 2 = (2-3A 2 )/(2A 2 ). This determines the critical 
value 

A c = (2.29) 

such that the edge vectors are timelike or null (incompatible inequalities) if |A| > A c , but 
spacelike (compatible inequalities) if |A| < A c . Note that the value of A that arises from 
dimensionally reducing 5-dimensional vacuum gravity down to 4 dimensions is A = and 
exceeds the critical value. This "explains" the conclusion reached in |14[] that the gravity- 
dilaton-Maxwell system obtained by Kaluza-Klein reduction of 5-dimensional gravity is 
oscillatory. 

We shall assume from now on that |A| < A c and that the Kasner exponents fulfill the 
above inequalities. For later use, we choose a number o > so that, for all x G Uq, 
o < 2pi — \A, a < 2pi + A A and a < pi/2. Reduce e if necessary so that e < a/7. If e 
is reduced, it may be necessary to shrink Uq so that the conditions imposed in section |27^ 
remain satisfied. In section [2.5| it is assumed that e and Uq are such that the conditions 
imposed in section [2.2| and the conditions imposed in this paragraph are all satisfied. 



2.5 Fuchsian system which is equivalent to the Einstein-matter 
evolution equations 

2.5.1 Rewriting of equations 

Theorem 3 in ]15[ (theorem 4.2 in preprint version), on which we rely for our result, states 
that a Fuchsian equation (i.e., as we mentioned above, an equation of the form ( |1.5D where 
A satisfies a positivity condition and / is regular, which includes a boundedness property) 
has a unique solution u that vanishes as t { 0, and furthermore spatial derivatives of u 
of any order vanish as t J. 0, as shown in Our goal is to recast the Einstein-matter 
evolution equations as a Fuchsian equation for the deviations from the velocity-dominated 
solutions. Thus, we denote the unknown vector u as 

u = K a b, if), u a , X , C, <Pab) (2.30) 

where the variables 7% etc. are related to the Einstein-matter variables by 



dab 


= 9ab + 9act b l b , 


(2.31) 


e c (l a b ) 


= t ^\ a ba 


(2.32) 


k a b 


= 9ac(°k C b + t- 1+ ^K C b), 


(2.33) 







(2.34) 




= t~ c u a , 


(2.35) 


td t ip + Pip 


= X, 


(2.36) 


S a 


= °S a + t^ a , 


(2.37) 


F a b 


= °F ab + tP<p ab . 


(2.38) 
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In the first of these equations t aCb is not the exponentiation of a matrix with components 
a c b such as occurs in (|2.13| ). The expression t aCb is for each fixed value of c and b the number 



which is t raised to the power given by the number a c b (defined below). In equations (|2.31| ) 
and ( p.33| ) there is no summation on the index b (but there is a summation on c). In 



( [2.38] ) (fab is a totally antisymmetric spatial tensor, which contributes three independent 
components to u. This assumption is consistent with the form of the evolution equation 
for ip ab , equation ( p.46|) below. The exponents appearing in (|2.31|) - (|2.38|) are as follows. 



Define ao — 4e, (3 — e/100 and £ = e/200 (where e is the same (small) quantity which 
entered the definition of the auxiliary exponents q a in section |272| and which was further 
restricted at the end of section |2T4l ). All of these quantities are independent of t and x. 
Finally define 

a a b = 2 m&x(q b - q a , 0) + a = 2g max{ajf)} - 2q a + a . 

Note that the numbers a a b are all strictly positive. In the second definition of a a b we 
have used the fact that the g a 's are ordered. The role of a a b is to shift the spectrum of 



the Fuchsian-system matrix A, in equation ( |1.5| ), to be positive. It is not clear to what 



extent the choice of a a b is fixed by the requirement of getting a Fuchsian system. It seems 
that the (triangle-like) inequality (42) of [|I5| (inequality (5.9) in preprint version) is a key 



property of these coefficients. We shall further comment below on the specific choice of 
a a b and its link with the BKL-type approach to the cosmological behaviour near t = 0. 

When writing the first-order evolution system for u we momentarily abandon the re- 
striction that g ab and k ab be symmetric, as in iJlSj . Thus we need to define g ab , and we 
do so by requiring that g a bg bc = 5 a c - This implies that g ab g bc = 5 a c . We lower indices on 
tensors by contraction with the second index of g ab and also raise indices on tensors by 
contraction with the second index of g ab . This choice is so that raising and then lowering 
a given index results in the original tensor, and the same for lowering and then raising an 
index. The position of the indices on quantities appearing in u and other such quantities 
is fixed. Repeated indices on these quantities imply a summation. On the other hand, as 
we already mentioned above, one qualifies the summation convention by insisting that in- 
dices repeated only because of their occurrence on p a , q a , a a b and other such non-tensorial 
quantities should be ignored when determining sums. 

Substituting ( |2.31| ) - ( |2.38| ) in the evolution equations yields equations of motion for u 
of the form (II. 51) 



td t r b + a\r b + 2k\ - 2{t°k\) 1 c b + 2r c {t Q k\) = -2t aa c +ac »- aa »r c K c b (2.39) 

t d t X a bc = fe c (t d tl \) + C t f e c ( 7 a J (2.40) 

td t K a b + a a b K a b - (t°k a b )(tiK) = t ao (tTK)K a b + t 2 - aab ( s R a b - M a b ) (2.41) 

td t tp + Pip-x = (2.42) 

td t oo a = t<{e a ( X ) + (C - (3)eaW} (2.43) 
td t x + Px = t a °- (3 (tTK)(A + t (i x)+t 2 -' 3 g ab V a V b (p + t 2 - (: V a u a 

+ t2 ^{^9ab£ a £ b e- x <f ) - ^g ab g ch F ac F bh e^} (2.44) 
td t e + PC = t^{e b {J-gg ac g hh F ch e^) 

Hfl9 ac + \m c )^9 hh F ch e^} (2.45) 
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t d^ab + Pfab = t 1 - /3 {-2e [a (^ ]c ^ c e-^) + fcJ-g^e-*} (2.46) 



All the quantities entering these equations have been defined, except s R a b - This is done 



by taking the Ricci tensor of the symmetric part of g ab jl5| . More explicitly, R a b 
g ac s R cb , with 

S d _ 4-qa+qb S p h _ 4.q a +qb SB ( s V h \ _ B ( s V h \ _i_ s V i s V h _ s V i s V h _i_ ■ S F' 1 1 

J^ab — «- -fi-ahb — 6 l e /A 1 afeJ e al 1 ~r 1 aft 1 /ii 1 hft L ai ^~ J ah L ibh 



(2.47) 



and the connection coefficients in the frame {e a }, 

X 1 

= - s g ch {e a (9(bh)) + e b (g {ha) ) - e h (g {ab) ) - g (ia) f bh - ~g {b i)~f ah } + -f ab - (2.48) 



Here. 

Sgab ig 

defined as the inverse of g( a b)- Once it is shown that the tensor g ab in equation 
( p.31| ) is symmetric, then it follows that s R a b = R a b and that equations (|2.39|) - ( p.46|) are 
equivalent to the Einstein-matter equations. 



2.5.2 The system ( ^739| )-( p^6|) is Fuchsian 



A good deal of the work needed to show that equation ( |1.5|) (as written out in equations 
( p.39| ) - (|2.46|) ) is Fuchsian was done in [fTql , in the massless scalar field case considered 
there. The form of the velocity-dominated evolution and the form of the function u are 
the same in the two cases except for the crucial addition of new source terms and new 
evolution equations involving the Maxwell field. The presence of the new components does 
not alter already existing parts of the matrix A, nor already existing terms in /. The 
difference between A here and A in the massless scalar field case considered in jnj is that 
here there are additional rows and columns, such that the only non-vanishing new entries 
are on the diagonal and strictly positive. Therefore the argument in [|Tj| that their A 
satisfies the appropriate positivity condition implies that our A satisfies the appropriate 
positivity condition. 

On the other hand, it is crucial to control in detail the new source terms in /, connected 
to the Maxwell field, which were absent in ||15|| . It is for the study of these terms that the 
results of [|IJ, and in particular the inequalities ( |2.24j ) which were shown there to guarantee 
that Maxwell source terms become asymptotically subdominant near the singularity, be- 
come important. Recall that the crucial criterion for the source f(t,x,u,u x ) is that it be 
0(t s ), for some strictly positive 5. In regard to this estimate, we use the notation "big 0," 
and "small o" as follows. Given two functions F(t,x,u,u x ) and G(t,x,u,u x ) we use 
the notation F ^ G, to denote that, for every compact set K, there exists a constant C 
and a number t > such that \F(t,x,u,u x )\ < C\G(t, x, u, u x )\ when (x, u, u x ) G K and 
< t < to (see Definition 1 in (|Tj|). If G is a function only of t (e.g. a power of t), then 
we replace F ^ G with F = 0(G). If f(t,x,u,u x ) = 0(t 5 ), then by reducing the value of 
6 (keeping it positive) we have that f(t,x,u,u x ) = o(t s ) with a "small o" which denotes 
that / /t s tends to zero uniformly on compact sets K as t — » 0. 

The new source terms involving the Maxwell field are: the last four terms in M a b (see 
equation ( [2.1|) ), the last two terms on the right hand side of equation ( |2.44| ) and the terms 
of the right hand sides of equations ( 2.45 ) and ( 2.46 ). 
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The calculation of the estimates starts in the frame, {e a }, defined in section p.2| . For 



more details concerning the basic estimates, we refer the reader to JT^]. In the frame {e a } 
the Kasner-like metric is (cf. (|2.13| )) 



°~g ab 



0~ab 



7. 



9ac(t 



t \ 2(K—Q) 
t ^-2(K-QP a 

T 



V(*o) 



(2.49) 
(2.50) 



where the matrix Q is the diagonal matrix Q a b = q a S a b which commutes with K. With our 
choice of frame, °g a b(to) = 5 a b and °g ab (to) = 5 ab . In Lemma 2 in |15|] (lemma 5.1 in preprint 
version), the form of ( |2.49| ) and ( |2.50| ) is considered and it is shown that °g ab = 0{t~ e ) and 
°g ab = 0(t~ e ). It is useful to write down expressions for the proposed metric and extrinsic 
curvature in the frame {e a }. The components in terms of this frame are 



9ab 
k a b 



°9ab+ °~g ac t &c »7 c b , 

~9ac(°k C b + t- 1+&C b K C b ). 



Here, a a b 



oc a b + q a — qb = \q a — Qb \ +«o is symmetric in a, b, a a b = a a . To get an estimate 
for the inverse metric, we note first that the inverse of g a c°9 cb is given by g ca0 g c b- Thus it 
is possible to express the latter quantity algebraically in terms of °g ab and r y a b . Now define 



Tb = 

which, from what we just observed, 
quantities and 7%. Then one has 



-t b [6 b 



~ac0~ \ 

9 9cb), 



(2.51) 



can be expressed algebraically in terms of known 



~ab 



_0 ~ab + t c 



^ a c °g cb . 



(2.52) 



As a consequence of an argument given in |15] which uses the (triangle-like) inequality (42) 
of that paper ((5.9) in preprint version) and the matrix identity preceding it, this exhibits 
7 a 6 as a regular function of 7%. In particular, if it is known that 7% is o(l) then the same 
is true of 7%. 

To better grasp the usefulness of the introduction of the exponents a a b and a a b , and 
the link of the Fuchsian estimates with the approximate estimates used in the BKL-like 
works, let us consider more closely the simple case where all the Kasner exponents are 
distinct (Case III). In this case p a = q a and one can diagonalize the Kasner-metric, so 
that, in the rescaled frame e a , we have simply (for all t < t ) °g a b(t) = 8ab- In such a case, 
the BKL-type estimates would be obtained (in the time-dependent rescaled frame e a ) by 



approximating the exact metric by its Kasner limit, i.e. 



~BKL 



simply g ab 



it) = 5 a 



>,. By contrast, 

the estimates of the Fuchsian analysis are made with the exact metric, g a b(t) = 5 ab +t aa b^ a b) 
in which / y a b , being part of u, is considered to be in a compact set and hence is bounded. 
As the diagonal a a a = > 0, we see that (in the frame e a ) the diagonal components of 
the "Fuchsian" metric asymptote those of the "BKL" metric, and that both are close to 
one. Concerning the non-diagonal components (in the frame e a ) of the "Fuchsian" metric 
we see that they are constrained, by construction (i.e. by the choice a a b = \q a — q b \ + a ), to 
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tend to zero faster than p 9a ~ 9b K This closeness between the metrics used in the two types 
of estimates explains the parallelism between the rigorous results derived here and the 
heuristic estimates used in BKL-type works. If we come back to the general case where 
the Kasner metric cannot be diagonalized in an analytic fashion, the optimal estimates 
become worse by a negative power of t (coming from the estimate of the matrix difference 
2(K—Q) in equations ( |2.49| ), (|2.50 ) above). The proposed metric in the frame {e a } satisfies 



then 

9ab =< t^-^ and g ab ■< t^-^H. 
The proposed inverse metric in the adapted frame is 

g ab = O g ab + t a« c -a c O g cb 

The proposed metric in the adapted frame satisfies 

9ab =< t 2q ^ a ^ e and g ab ± r 29mi ^»~ e . (2.53) 

Estimates of spatial derivatives of the proposed metric are also needed. 

e c (9ab) =< and e c (g ab ) ^ 

e c (9ab) =< t 2 *—^.*}-*" 6 and e c (# afe ) r< fr- 2 *"**^ - ' -6 (2.54) 

for some strictly positive 5. 

The determinant of the proposed metric also appears in some of the new source terms. 
From (|2.13|) , the form of °g a b{to) and trff = 1, one gets °g = t 2 . From ( 2.49|) and 



°9ab(to) — $ab one gets °g = 1. The expression for the determinant is a sum of terms of 
the form g a b9cd9ef, such that in each term, each index, 1, 2, 3, occurs exactly twice. From 
the Kasner relation for the g a 's and the relation between the two frames, it follows that 
g = t 2 g. Considering the form of the various expressions, one then obtains 1/g — 0(t~ 2 ), 
^g = 0(t), 1/y/g = 0(t- l ) } and 1/Jg-l/^j = 0(r 1+ ^) = 0(r 1+e ). Spatial 
derivatives of the determinant also appear in /. Considering the form of g — °g and that 
e a (°g) = 0, it follows that e a (g) = 0(t a °- s - 3€ ), and 

e a (g)=0(t 2+a °- s -^). (2.55) 

Finally, 

e a (g- 1/2 ) = -g^ = 0(r^^). 

Let us now consider the new source terms in /, beginning with the last four terms of 
t 2 ~ aa ' b M a i ) . To estimate the contributions of £ a and F ab it is sufficient to note from ( |2.37| ) 
and ( p.38| ) that £ a = 0(1) and F ab = 0(1). Then we get, using the definition of a a b and 

t 2 - a ^-{g hc E a E c - U a b g ch S c S h }e-^ 
9 2 

J- 2 W{a,6}+ 2 9"+ 2 9max{b,c}-^- a 0-e _|_ £+2<2 max { c>h } - XA-a -e 
c c,h 

19 



^ t 2qx-XA-a -e = q^- 



0(t 6 



t 2 - aa "{g ac g M F ch F bl - ~6 a b g ch g ij F ci F hj }e x * 

-<; ^2-2g max{aji) }+2(}a-2(j min { aiC }-2g mill { hii }+AA-ao-2e 

_|_ ^2-2g min { c j l y-2q rnin { ij y+\A-a -2£ 

c,h,i^c,jj^h 



0(r 



0(f 



for some strictly positive 5. The crucial inputs in getting these estimates are the inequalities 
Q2.24Q . We recall also that the quantity a (linked to (|2.24|) being satisfied) was introduced 



at the end of subsection |2.4| , The estimate of the last two terms on the right hand side of 
(EH) ^ 



t^-g ab S a S b e-^ 
9 

t 2 -?g ab g ch F ac F bh e x * 



o(r p - 2£+a ) 



o(t 5 ) 



The right hand side of fl2~45l) is 0{t a °-^~ 5 ~ be+CT ) = 0{t s ). The right hand side of (gig) 
is (9(£ Q o-/3-5-4e+o-^ — 0{t 5 ). The other terms which occur in / were estimated in | 15| , 
resulting in / = 0(t 5 ). 

To show that we indeed have a Fuchsian equation, we need to check not only that 
/ = 0(t s ), but also that d u f = 0(t 5 ) and d Ux f = 0(t s ), along with other regularity 
conditions JT5], 2~3|. In JTB| it is shown that / is regular with equation (31) in that paper and 
the remarks following equation (31). In our case there is a factor involving the determinant 
of the metric in various of the terms in / which are not present in the case considered in |T5| . 
The discussion surrounding equation (31) in |TJj applies to our case as well, even for terms 
in / containing g ±1//2 '. The Kasner-like contribution is the leading term, and this function 
of t and x can be factored out. What is left is of the form w(t, x, u, u x )(l + h(t, x, u))^ 1 ^ 2 , 
which is analytic in h at h — 0. The conditions listed following equation (31) hold. Thus 
we conclude that ( |1.5| ) as written out in ( |2.39| ) - ( 2.46 ) is a Fuchsian equation. 



2.5.3 Symmetry of metric 



It remains to show that g ab is symmetric, so that equation ( |1.5|) as written out in (|2.39|) 
- ( [2.46| ) is equivalent to the Einstein-matter evolution equations. The structure of the 
argument is the same in any dimension and so it will be written down for general The 
number of distinct eigenvalues of K a b is maximal almost everywhere. Thus it is enough to 
show that g\ ab -\ and k\ a u vanish in the case that the Kasner-like metric is diagonal, since 

6 The argument for the symmetry of the metric in fl5|| is not valid as written since some terms were 
omitted in the evolution equation for the antisymmetric part of the extrinsic curvature. The correct 
equation is 

d t (kab ~ ha) = (trk)(k ab - k ha ) - 2(k ac k c b - k bc k c a ). 

In the following a proof of the symmetry of the tensors g ab and k ab is supplied with the help of a different 
method. 
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then by analytic continuation they vanish on the entire domain. We therefore consider the 
case that the Kasner-like metric is diagonal. 

The redefinitions (|2.31|) , Q2.33Q from the variables g ab , k ab to the variables 7 a &, n a b were 



viewed in the previous subsections as a change between variables with no particular sym- 
metry properties in their indices (18 on each side). One can, however enforce g[ ab ] = by 
assuming that 7% is symmetric and vice-versa. Indeed, under our diagonality assumption 
for the Kasner-like metric, g a b{f) = 8 ab + t aab j a b where a a b = \q a — q b \ + a is symmetric in 
(a,b). Accordingly, imposing the symmetry 7% = 7 b a algebraically ensures the symmetry 
of g ab . Similarly, one can enforce k[ ab ] to vanish by imposing consistent constraints on 
inserting ( [2.3 1| ) into ( |2.33| ) (with the velocity-dominated solution diagonal) and writing 
out the constraint k ab — k ba = gives the following condition on K a b 

«°& - K b a - 7 > 6 + j b aPa + ( 7 > c b - 7 V a ) = 0, (2.56) 

with a (afe)c = 2p max { ajC } + 2p max{b>c} - 2p max{ayb} - 2p c + a . These conditions show that 
there are only six independent components among the K a b , which can be taken to be those 
with a < b. This is because, the relation ( |2.56|) can be solved uniquely for the components 



K a b with a > b, given the other ones, at least for t small. That this is true can be seen 
as follows. Rearrange the equations fl2.56j ) so that the terms containing K a b with a > b 
are on the left hand side and all other terms are on the right hand side. The result is 
an inhomogeneous linear system of the form A(t, x)v (t, x) = w(t,x) where A(t, x) and 
w (t, x) are known quantities and v denotes the components K a b with a > b which we want 
to determine. Furthermore A(t, x) — I + o(l), where I denotes the identity matrix. It 
follows that A(t) is invertible for t small, which is what we wanted to show. The solution 
n a b (a > b) remains moreover bounded when 7 b a and K a b are in a compact set. We shall 



assume from now on that ^ a b is symmetric and K a b constrained by ( |2.56|) , so that symmetry 



of the metric is automatic. The redefinitions (p ■ 3 1|) , (|2.33|) from g ab , k ab to ^ a b) K a b can 



now be viewed as an invertible change of variables, from 12 independent variables to 12 
independent variables. We can also clearly assume X a bc in (|2.32|) to be symmetric in a, b. 

With these conventions, there are less components in u than in the previous subsections. 
The independent components can be taken to be 7%, K a b and X a bc with a < b, together 
with the matter variables. An independent system of evolution equations is given by ( |2.39| ) 
- (|2.41|) with a < b for the gravitational variables, and the same evolution equations as 



before for the matter variables. These evolution equations are equivalent to all the original 
evolution equations, since the equations ( |2.39| ) - ( |2.41| ) with a > b are then automatically 
fulfilled, as can be shown using the fact that the Einstein tensor and the stress-energy tensor 
are symmetric for symmetric metrics. To see this it must be shown that given a symmetric 
tensor S ab , the vanishing of S a b = g ac S cb for a < b implies that S ab = 0. Consider the 
linear map which takes a symmetric tensor S ab , raises an index, and keeps the components 
of the result with a < b. This is a mapping between vector spaces of dimension d(d+ l)/2 
and can be shown to be an isomorphism by elementary linear algebra. This proves the 
desired result. 

Now, this reduced evolution system is also Fuchsian. This follows from the same reason- 
ing as above, which still holds because all components of u, including the non-independent 
ones, can still be assumed to be bounded. Therefore, there is a unique u that goes to 
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zero, which must be equal to the one considered in the previous subsections. The metric 
considered previously is thus indeed symmetric. 



2.5.4 Unique solution on a neighbourhood of the singularity 

Given an analytic solution to the velocity-dominated evolution equations on (0, oo) x 
E, such that inequalities ( |2.24| ) are satisfied, we now have a solution u to a Fuchsian 
equation (and a corresponding solution to the Einstein-matter evolution equations) in 
the intersection of a neighbourhood of the singularity with (0, oo) x Uq where Uq is a 
neighbourhood of an arbitrary point on S. These local solutions can be patched together 
to get a solution to the Einstein-matter evolution equations everywhere in space near the 
singularity. It may seem like there could be a problem patching together the solutions 
obtained on distinct neighborhoods with non-empty intersection because the Fuchsian 
equation is not the same for different allowed choices of e and adapted local frame. The 
construction is possible because different allowed choices of e and local frame result in a 
well-defined relationship between the different solutions u which are obtained, such that the 
corresponding Einstein-matter variables agree on the intersection (up to change of basis). 
It therefore follows that given an analytic solution to the velocity-dominated evolution 
equations on (0, oo) xE, such that inequalities (|2.24|) are satisfied, our construction uniquely 



determines a solution to the Einstein-matter evolution equations everywhere in space, near 
the singularity. 



2.6 Einstein-matter constraints 
2.6.1 Matter constraints 

The time derivative of the matter constraint quantities (the left hand side of equations ( |2.2|) 
and (|2.3| )) vanishes. If the velocity-dominated matter constraints are satisfied, the matter 
constraint quantities are o(l). A quantity which is constant in time and o(l) must vanish. 
Therefore the matter constraints are satisfied. 



2.6.2 Diagonal Kasner metrics 

It remains to show that the Hamiltonian and momentum constraints are satisfied, that C 
and C a , defined in (|1.9| ) and ( |1.10| ), vanish. Since we now have a metric, g^, it follows that 



V M G^ = 0. Since the matter evolution and constraint equations are satisfied, it follows 
that V ll T^ 1 ' = 0. From the vanishing of the right hand side of ( |1.12 ) and the vanishing of 



the covariant divergence of both the Einstein tensor and the stress-energy tensor, it follows 
that 

d t C = 2{trk)C -2V a C a (2.57) 
d t C a = {tik)C a - l -V a C. (2.58) 

Now define C = t 2 -^C and C a = t l ~^C a , with < r/ 2 < ^ < (3. 

tdtC + r^C = 2(1 + t tr k)C -2t 2 - Vl+m V a C a (2.59) 
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td t C a + V2 C a = (l + ttrk)C a --t^V a C 



(2.60) 



On the right hand side of Q2.59D and Q2.60| ) C and C a are to be considered as components 
of u — (C, C a ). If it is shown that ( |2.59| ) and ( |2.60| ) is a Fuchsian system, then there is a 
unique solution u such that u = o(l). It is clear that u = is a solution to (|2.59|) and Q2.60|) . 
If it is shown that C = o(l) and C a = o(l), (i.e. that C = o(t" 2+r ' 1 ) and C a = o(t" 1+r ' 2 )), 
then they must be this unique solution. Furthermore, it is sufficient to consider the case 
that the Kasner-like metric is diagonal, since the number of distinct eigenvalues of K a b is 
maximal on an open set of E. If the constraints vanish on an open set of their domain, 
then by analytic continuation they vanish everywhere on their domain. 

Therefore we consider the case that the Kasner-like metric is diagonal and show first 
that 



1 + ttrJfe 
V a C a 



0(t 5 ) 



(2.61) 
(2.62) 



(when C a is bounded) so that the system (]2.59|) , (|2.60| ) is Fuchsian (the complete regularity 
of f(t,x,u,u x ) defined by fl2.59Q and ( |2.60| ) can be easily verified); and second, that 



C 

c n 



o(t- l+r > 2 ) 



(2.63) 
(2.64) 



Some facts which will be used to show this follow. Consider indices a G {1,2, 3}. The 
following inequalities hold for some positive integer n and for real numbers, q a , ordered 
such that if a < b, then q a < q^. (In later sections we define ordered auxiliary exponents, 
q a , for a G {1, ■■■,d}, for arbitrary fixed d > 2. Then ( |2.65| ) - ( |2.67| ) hold more generally 
for indices in {1, d}.) 

n 

Vao + Yl kai-i - I - ^ > (2.65) 



i=l 



: 



maxjafcjaj} 



i=l 
n 



-Qao+Yl I -Qai\- Qan > -H 



minjaj.,^'} 



(2.66) 
(2.67) 



i=l 



The latter two inequalities hold for any k, j in {0, . . . , n}. 

In the case that the Kasner-like metric is diagonal, q a = p a . The metric in the frame 
{e a } is °g ab = Sab, 



9ab 

~ ab 



Pa-Pb\ 



5*b + t a\-a b ^ t \p a -Pb\_ 



5 ab + t«^ a b =< t 

a 

b 



The extrinsic curvature satisfies t °k a b = —5 a bPb, 



tk\ 

t (k a b — °k a b ) 



-5 a bp b + t a ^ a b , 
b- 
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and t tr °k = —1 



ttrfc = -l+t ao tTK, (2.68) 

t 2 {(tr£;) 2 -(tr £;) 2 } = 0(t a °). (2.69) 

The following estimates will also be useful. 

_ h a , b b -L b a ,°b b — n.-2+a a _ ±~2+a a b +a b a a b 

= O(r 2+ao ), (2.70) 

and 

e a (trk - tr°Jfe) = e a (t- 1+ao tr k) = 0(r 1+Q °). (2.71) 
The structure functions of the frame {e a } are 

fib = t p ^~ Pb f c ab -\nte a (p b )t~ Pa 5 c b + \nte b (p a )r Pb 6 c a 
It is convenient to have an estimate of T c ab , the connection coefficients (|2.48|) in the frame 



Term A: 


~9 ch e a (9bh) ^J2 tlPc ' 

h 
h 

~g ch e h (~gab)±Y, tlPc ~ 

h 

rhafl < E ^ 

h,ij^b 


~Ph\ 


-Pa+\Pb~Ph\ 


-5 

i 


(2.72) 


Term B: 


~Ph\ 


-Pb+\Ph~Pa\ 


-5 

> 


(2.73) 


Term C: 


~Ph\ 


-Ph+\Pa~Pb\ 


-S 


(2.74) 


Term D: 


-Pi|-+ 


\Ph-Pa\+Ph~ 


-Pb-Pi-S 

? 


(2.75) 


Term E: 


~g ci g bh f* l E tlPc ~ 


-pM 


\Pb~Ph\+Ph~ 


-Pa-Pi-S 

) 


(2.76) 


Term F: 










(2.77) 



The difference between the connection coefficients for the metric g ab and those for the 
Kasner-like metric is Ar^ b = — °r^ 6 . It is useful to have the estimates 

f I = l~g ab e c (g ab ) + ] a ac < t- p <- 5 , 

and 

Alt = l? b e c (g ab ) =< t- p ^~ s . (2.78) 
2.6.3 Momentum and Hamiltonian constraints 

First, we show (|2.61|) and (|2.62|) . From equation ( [2.68D , 1 +ttrk — 0(t a °). Similarly, we 
can estimate V a C a , 

V°C a = g ab V a C b 

= §*{r p °e a (C b t-*>)-T c ab C e t-*>} 
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The first term is 

g ab f-Pa. e a (Cbt~ Pb ) -< t\P a ~Pb\-Pa-Pi-& ^--2p min { a ,6}-5 

From (|2.72| ) - ( |2.76|) the second term is 

g ab t c ab c c r Pc < r 2ps - s . 



[2.79) 



(2.80) 



From (|JD and (ET80D , the desired estimate, V a C a = 0{t- 2+ ^ 2 ) is obtained. Thus, the 
system ( |2.59| ), ( |2.60| ) is Fuchsian. 

Next we turn to (|2T63|) and |T64T ). A term that appears in the momentum constrain; 
is V a fcV The estimate is needed in the adapted frame, and the covariant derivative is 
calculated in the frame {e a }. This adds a factor of t Pb , 

V a /c a fe = {e a (k a b ) + T a ac k c b — T c ab k a c } t Pb 

Furthermore, the quantity whose estimate will be required is the difference between this 
term and the corresponding term in the velocity-dominated constraint, 



b 



°V n °k a , 



{e»(ky - »k\) + ATI + T a ac (k c b - »k c b )}t 

F c b a +Pb _L °F C °b a +Pb 
1 _i.fi/ r L -|- 1 r , h ft, „ L 



(2.81) 
(2.82) 



The right hand side of (|2.81|) is O(t~ 1+ao ~ s ). The terms in line (|2.82| ) originating from 
Term E of the connection coefficients (see ( j2.76| )) are cancelled in the sum, due to the 
antisymmetry of and the symmetry of k at and °k at . For estimating the rest of the 
terms in line ( |2.82|) , it is convenient to rewrite this line as, 



r c h a +Pb _l °r c °b a +pi> 



AT c ab »k a c t Pb -T c ab (k a c -»k a c )t Pb 



with 



Ar 



ab 



h 

i h hi 



(2.83) 

(2.84) 
(2.85) 

(2.86) 

(2.87) 



hi 



The terms in line ( p.87 ) need not be considered since they originate from Term E of the 
connection coefficients and as stated above the contribution from this term is cancelled by 
terms in A = Y c ab ®k a c t Pb . So considering only lines (|2.84|) - fl2.86|) , the first term on the 
right hand side of (|2.83|) is 



AT a ab p a r 1+Pb = O(r 1+ao ~ 5 ) + terms which are cancelled by A. (2.88) 

Since the terms in the sum come with different weights, p a , ( [2.78|) cannot be used in ( |2.88| ). 
But the estimate is straightforward. For example, the term in ( 2.88 ) originating from 
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the 3rd term in line ( ggg) ^ Y, a , h ,i t ~ lMpa ~ Ph ^ Pi ~ Pa ^ i ~ Ph+2a °~ S = 0(r 1+Q °- 5 ). Finally 
consider the rest of the right hand side of (|2.83|) , 



~9 ch {e a (~g bh ) + e b (9ha) ~ e h (g ab ) - ~gJl h - ~gj l ah } + f c ab \ t &a ^ a c r 1+p \ (2.89) 



For all terms except the 5th term in ( |2.89| ), the estimate, O(t~ 1+ao ~ 5 ) can be obtained 
from ( |2.72p - ( |2.77| ). The fifth term originates from Term E of the connection coefficients, 
and was already considered above. Therefore 



V a fc a 6- °V a °fc a - 



0(f 



-l+a -8\ 



(2.90) 



Next the matter terms are estimated. For the Hamiltonian constraint, an estimate of 
p — °p is needed. 

(^0) 2 -(^V) 2 = {2d t G <p + r 1+(3 {(3^ + td^)}r l+f3 {(3^ + td t ^) 

= {2A + ^{(3tfj + td t ^)}r 2+f3 {^ + td^) 

= (r 2+ ^), 
g ab e a (4>)e b (<f>) ± r 2p °- s - e = o(r 2+r *), 



-g ab S a S b e- x ^ 

g 

g ab g ch F ac F bh e x ^ 



0(r 2 - e+a ) 
o(r 2 ~ 2e+a ) 



o(r 2 ^). 

-- o{t- 2+r ^ 



Therefore, 



p- °p = o(t- 2+ ^). (2.91) 
The difference between the matter terms in the momentum constraint and in °C a is 

: 0(r 1+ e- s ), 



-^e a (0) + 9iVe a (V) = -d t °<Pe a {t^)-d t {t^)e a {<P) 

( ] _J ^0p60 p , 1 

{—p 7fj= t £ab + — 

\[g v°9 Vd 



°S b0 F ab + ^(S b F ab -°S b0 F ab ) = o(r 1+ " 2 ). (2.92) 



Estimates related to the determinant which are relevant to ( |2.92| ) can be found immediately 
preceding equation (|2.55|) . From the estimates just obtained, 



o(t 



(2.93) 



From R = 0(t- 2+a °) (shown in pf) and from °C = 0, (|2T70| ), ( ggg) , ( ggg) and the 
relative magnitude of the various exponents, it follows that C = o(t~ 2+Vl ). From °C a = 0, 
( |2.90| ), (|2.71| ), ( p.93| ) and the relative magnitude of the various exponents, it follows that 

c a = (t- i+ n 

Since ( |2.63| ) - (|2.64 ) are satisfied, the Hamiltonian and momentum constraints are 
satisfied. 
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2.7 Counting the number of arbitrary functions 



The number of degrees of freedom of the Einstein-Maxwell-dilaton system in 4 spacetime 
dimensions is 5: 2 for the gravitational field, 2 for the electromagnetic field and 1 for the 
dilaton. Hence, a general solution to the equations of motion should contain 10 freely 
adjustable, physically relevant, functions of space (each degree of freedom needs two initial 
data, q and q). This is exactly the number that appears in the above Kasner-like solutions. 

• The metric carries four, physically relevant, distinct functions of space. This is the 
standard calculation 

• The scalar field carries two functions of space, A and B. 

• The electromagnetic field carries six functions of space, °S a and °F a b. These are 
physically relevant because they are gauge invariant, but they are subject to two 
constraints, leaving four independent functions. 

A different way to arrive at the same conclusions is to observe that the respective number 
of fields, dynamical equations and (first class) constraints are the same for the velocity- 
dominated system and the full system. Hence, a general solution of the velocity-dominated 
system (in the sense of function counting) will contain the same number of physically 
distinct, arbitrary functions as a general solution of the full system. This general argument 
applies to all systems considered below and hence will not be repeated. 



In |L5| a different way of assessing the generality of the solutions constructed was used. 
This involved exhibiting a correspondence between solutions of the velocity-dominated 
constraints and solutions of the full constraints using the conformal method. That method 
starts with certain free data and shows the existence of a unique solution of the constraints 
corresponding to each set of free data. It is a standard method for exploring the solution 



space of the full Einstein constraints |42[ and in |1J| it was shown how to modify it to 
apply to the velocity-dominated constraints. While it is likely that the conformal method 
can be applied in some way to all the matter models considered in this paper, the details 
will only be worked out in two cases which suffice to illustrate the main aspects of the 
procedure. These are the Einstein-Maxwell-dilaton system with D = 4 (this section) and 
the Einstein vacuum equations with arbitrary D > 4 (next section). Even in those cases no 
attempt will be made to give an exhaustive treatment of all issues arising. It will, however, 
be shown that the strategies presented for solving the velocity-dominated constraints are 
successful in some important situations. 



The procedure presented in the following is slightly different from that used in [15 



Even for the case of the Einstein-scalar field system with D — 4 it gives results which 



are in principle stronger than those in [15| since they are not confined to solutions which 
are close to isotropic ones. In the presence of exponential dilaton couplings a change of 
method seems unavoidable. One part of the conformal method concerns the construction 
of symmetric second rank tensors which are traceless and have prescribed divergence from 
the truly free data. In this step there is no difference between the full constraints and the 
velocity-dominated ones. An account of the methods applied to the full constraints in the 



case D = 4 can be found in [[42]. (These arguments generalize in a straightforward way 



to other D. It is merely necessary to find the correct conformal rescalings. For D > 4 
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and vacuum these are written down in the next section.) In view of this we say, with a 
slight abuse of terminology, that the free data consists of a collection g ab , k ab , H, <fi, <fi t , £ a , 
F ab where g ab is a Riemannian metric, k ab is a symmetric tensor with vanishing trace and 
prescribed divergence with respect to g ab , H is a non-zero constant, <fi and <j> t are scalar 
functions and S a and F ab are objects of the same kind as elsewhere in this section. All 
these objects are defined on a three-dimensional manifold. Next we introduce a positive 
real- valued function ui which is used to construct solutions of the constraints from the free 
data. Define g ab = w 4 £u, Kb = uj~ 2 k ab + Hg ab , <p t = u~ 6 4> t .JThe objects g ab , k ab , <p, <p t , 
£ a and F ab satisfy the constraints provided the divergence of k ab is prescribed as u 6 j a and 
u> satisfies a nonlinear equation which in the case of the full Einstein equations is known 
as the Lichnerowicz equation. In the case of the velocity-dominated constraints it is an 
algebraic equation. The Lichnerowicz equation is of the form 

A § uj - -Rgoo + J2 a i^ ai ~ \ h2 ^ = ( 2 - 94 ) 

8 i=l 4 

Here ot\ = 1, a 2 = — 3 and a 3 = —7. The functions a, depend on the free data and their 
exact form is unimportant. All that is of interest are that each a, is non-negative and that 
at any point of space a x = iff V a = , a 2 = iff the electromagnetic data vanish and 
a 3 = iff <p t and k ab vanish. Next consider the velocity-dominated constraints for d = 3. 
The analogue of the elliptic equation ( |2.94| ) is the algebraic equation 

bu- 7 - -H 2 u 5 = (2.95) 

Here b is a non-negative function which vanishes at a point of space iff (fit and k ab vanish. 
This can be solved trivially for uj > provided b does not vanish at any point since the 
mean curvature H is non-zero. For each choice of free data satisfying this non-vanishing 
condition there is a unique solution to of Q2.95p . 



In order to compare the sets of solutions of the full and velocity-dominated constraints 
in these two cases it remains to investigate the solvability of the elliptic equation (|2.94|) for 



uj. A discussion of this type of problem in any dimension can be found in ||43|| . We would 
like to show that for suitable metrics on a compact manifold the equation for uj always 
has a unique solution, i.e., the situation is exactly as in the case of the velocity-dominated 
equations. The problem can be simplified by the use of the Yamabe theorem, which says 
that any metric can be conformally transformed to a metric of constant scalar curvature 
—1, or 1. In the following only the cases of negative and vanishing scalar curvature of 
the metric supplied by the Yamabe theorem will be considered. A key role in the existence 



and uniqueness theorems for equation ( |2.94|) is played by the positive zeros of the algebraic 



expressions x + 8 J2i=i a>iX ai — 6H 2 x 5 and 8 J2i=i UiX ai — 6H 2 x 5 . Provided J2 a =i a i does not 
vanish anywhere it is possible to show that each of the algebraic expressions has a unique 
positive zero for each value of the parameters. 

The significance of the information which has been obtained concerning the zeros of 
certain algebraic expressions is that it guarantees the existence of a positive solution of the 
corresponding elliptic equations for any set of free data satisfying the inequalities already 
stated using the method of sub- and supersolutions (cf. [43]). In the case of equation 
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( |2.94| ) uniqueness also holds. For in that case the equation has a form considered in for 
which uniqueness is demonstrated in that paper. The advantage of the three dimensional 
case is that there the problem reduces to the analysis of the roots of a cubic equation, 
a relatively simple task compared to the analysis of the zeros of the more complicated 
algebraic expressions occurring in higher dimensions. 

For the purpose of analysing Kasner-like (monotone) singularities it is not enough to 
know about producing just any solutions of the constraints. What we have shown is that (i) 
if the Kasner constraints are satisfied at time to, then they are propagated at all times by 
the velocity-dominated evolution equations; and (ii) if the Kasner constraints are satisfied, 
the exact constraints are also satisfied. It is also necessary to verify, however, that one can 
satisfy simultaneously the Kasner constraints and the inequalities necessary for applying 
the Fuchsian arguments, i.e., we must make sure that we can produce a sufficiently large 
class of solutions which satisfy the inequalities necessary to make them consistent with 
Kasner behaviour. Because of the indirect nature of the way of solving the momentum 
constraint (which has not been explained here) it is not easy to control the generalized 
Kasner exponents of the resulting spacetime. There is however, one practical possibility. 
Choose a spatially homogeneous solution with Abelian isometry group (for d = 3 this 
means Bianchi type I) which satisfies the necessary inequalities. Take the free data from 
that solution and deform it slightly. Then the generalized Kasner exponents of the final 
solution of the velocity-dominated equations will also only be changed slightly. If the 
homogeneous solution is defined on the torus T 3 then it is known that any other metric 
of constant scalar curvature has non-positive scalar curvature. Therefore we are in the 
case for which existence and uniqueness is discussed above. We could also start with a 
negatively curved Friedmann model. 



3 Vacuum solutions with D > 11 

The second class of solutions we consider is governed by the action ( |1.2|) , with D > 11. The 
d + 1 decomposition is as in section |1.6| , with the matter terms vanishing. The Kasner-like 
evolution equations are (|2.7|) and ( |2.8|) . The general analytic solution of these equations is 
given by ( |2. 12| ) and ( [2.13j ). To obtain this form, we again adapt a global time coordinate 



such that the singularity is at t — 0. We label the eigenvalues of K, p±, . . . ,Pd, such that 
Pa < Pb if o- < b. The eigenvalues again satisfy J2i=iPi — 1 ; coming from tr K = 1. As in 
the D = 4 case, in order to preserve analyticity even near the points where some of the 
eigenvalues coincide, while retaining control of the solution in terms of the eigenvalues, we 
introduce a special construction involving auxiliary exponents and an adapted frame. 

In higher dimensions, there are more possibilities to take care of, but the idea is the 
same as in the D = 4 case. Consider an arbitrary point xq G S. Let n be the number of 
distinct eigenvalues of K at x . Let be the multiplicity of i £ {1, . . . , n}, with p^ i 
such that pi, is strictly less than p^. if b < A4. Thus pa ■ ■ ■ ,PAi+mi-i are equal at x . For 
each integer a £ {A i: . . . , A4 + rrii — 1}, define 

Ai+rrii-l 
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on a neighbourhood of Xq, Uq. Note that if m; = 1, then qj^ = Pa v Shrinking Uq if 
necessary, choose e > such that for x G U , for a G {Ai, . . . , A{ + rrii — 1} and for 
b G {Aj, ...,Aj + m,j- 1}, if i = j, then \p a - p 6 | < e/2, while if i ^ j, \p a - p b \ > e/2. 

The adapted frame {e a } is again required to be analytic and such that the related 
frame {e a } is orthonormal with respect to the Kasner-like metric at some time to > 0, 
with e a = t~ Qa e a . In addition, it is required that e^, . . . , e^+mi-i span the eigenspace 
of K corresponding to the eigenvalues p^., . . . ,pA i+mi -i. Note that if rrii — 1 then e^. is 
an eigenvector of K corresponding to the eigenvalue q^.. The auxiliary exponents, q a , are 
analytic, satisfy the Kasner relation (J2 Qa = 1), are ordered (q a < q b for a < b), and satisfy 
q% > Pi, qd < Pd and max a \q a - p a \ < e/2. If q a ^ q b , then °g ab , °g ab , °g ab and °g ab all 
vanish, and the same is true with g replaced by k. 

The velocity-dominated constraints corresponding to the Hamiltonian and momentum 
constraints are °C = and °C a = 0, with °C and °C a as in equations ( 2.18|) and ( 2.19 ), 



with the matter terms vanishing. For the solution fl2.12p - ( |2.13| ) the velocity-dominated 



Hamiltonian constraint equation is equivalent to J^Pa = 1- Equations (|2.21|) - (|2.22|) are 
satisfied, so if the velocity-dominated constraints are satisfied at some t , then they are 
satisfied for all t > 0. 

For this class of solutions, the inequality [13|], 2p 1 +p 2 -\ +Pd-2 > 0, or equivalently, 



1+Pi -Pd-Pd-i > 0, (3.1) 



defines the set V which was referred to in the introduction. As shown in ||13|| , this inequality 
can be realized when the spacetime dimension D is equal to or greater than 11. As in our 
Maxwell archetypal example above, we expect that this inequality will be crucial to control 
the effect of the source terms (here linked to the spatial curvature) near the singularity. It is 
again convenient to introduce a number a > so that, for all x G U , 4cr < l+pi—pd—Pd-i- 
Reduce e if necessary so that e < a/{2d+ 1). If e is reduced, it may be necessary to shrink 
Uq so that the conditions imposed above remain satisfied. It is assumed that e and Uq are 
such that the conditions imposed above and the condition imposed in this paragraph are 
all satisfied. 

We again recast the evolution equations in the form ( |1.5|) and show, for D > 11, that 
( |1.5| ) is Fuchsian and equivalent to the vacuum Einstein equation, with quantities u, A 



and / as follows. Let u = (7%, X c e f, ^ h i) be related to the Einstein variables by (|2.31| ) 



- ( |2.33| ). For general d define a — (d + l)e and define a a b in terms of a as in section 
0. Let A and /, be given by equations ( |2.39| ) - ( |2.41| ), with M a b = 0. The argument 



that A in equation ( |L5|) satisfies the appropriate positivity condition is analogous to the 



part of the argument concerning the submatrix of A corresponding to (7, k) in 1 15]. A 



transformation to a frame in which °g ab is diagonal induces a similarity transformation of 



A. The eigenvalues of the submatrix are calculated in this representation in and the 
generalization of the calculation to integer d > 2 is straightforward. 

To obtain / = 0(t s ) requires the estimate t 2-0 " 6 s R a b = 0(t 5 ). The strategy used here 
is different from that used to estimate the curvature in (15| . The general problem is one of 
organization. There are many terms to be estimated, each of which on its own is not too 
difficult to handle. The difficulty is to maintain an overview of the different terms. The 



procedure in [15] made essential use of the fact that the indices only take three distinct 
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values and in the case of higher dimensions, where this simplification is not available, an 
alternative approach had to be developed. 

First s R a b is estimated by considering each of the five terms in the expression ( |2.47 ) 



These five terms are expanded by considering each of the six terms in fl2.48|) if the indices 



on s T^ b are distinct, but carrying out the summation before estimating s T° lb . There are 
thus 55 terms to estimate. While many of these terms are actually identical up to numerical 
factors, the ease with which each term can be estimated, using the inequalities (|2.65| ) - 
( [2.671 ), led to estimation of all 55 terms, rather than first combining terms. We do however, 
take into account that /- fc = if j = k for obtaining the estimates. 

Once an equation such as ( |1.5|) is shown to be Fuchsian, then it follows that spatial 
derivatives of u of any order are o(l). At the stage of the argument we are at here, we 
cannot assume u xx = 0(1). This means that t~^\ a bc must be used for eb(7 a c ) in places 
where a spatial derivative of eb(7 a c ) occurs. This makes a slight difference, compared to 
section in what estimate of the terms in the connection coefficients is used for the 
first and second terms of ( |2.47|) (t~ 5 is replaced by t - ^). There are additional differences 
from ( |2.72| ) - ( |2.77p , because there it is assumed that the Kasner-like metric is diagonal. 
The estimates °g a b = 0(t~ e ) and °g ab = 0(t~ e ), obtained in Lemma 2 of | 15| , hold in the 



case we are considering, so that ^ tl<?"-<?f>h £ and (see [[HJ) s g ab ^ t^ qa ~ qb ^~ e . This adds 
factors of t~ e to the estimate of terms in the connection coefficients. 
With these considerations, from (|2.48| ), 



S Kc = \ S ~9 ab e c (~ 9{ab) ) + tc 1 (3.2) 

Here we do not write out the estimates of all 55 terms, but instead give some examples, 
with a number designating which term of (|2.47|) is being considered (1 - 5), and a letter 



designating which term of (|2.48|) is being considered (A - F). Thus, for example, term 1C 



is 



t~ qa+Qb g ac eh( S g hi ei(g( c b))) ~^ ^ f-Qa+qb+\9a-qc\-qh+\ih-Qi\-Qi+\ic-qb\-S-Se-C 



c,h,; 



_^ £-2<ja+2<2 max{a b }-2g d -<5-3e-C ^ ^ 



Term 3E is 



^-Qa+qbgic S '~ik~ Jj^S^h ^-9o+gb+|g -gc|+|g»-gfc|+|g6-9j|+g>— 9c-9fe-9i-*-5e 

c,i,j,kj^c 

-<; ^-2g min{aiC }+2g f) -2g min { dift }-5-5e ^ 

In term 4, t~ qa+qb g ac s T^ b s T l ch , the terms resulting from expanding s T^ b are designated by 
small letters a - f , and those from T l ch are designated by capital letters A - F. Term 4dA 
is 

r*^V%i)/4V%(0<A*)) 

£-qa+qb+\qa-qc\ + \qh-qi\ + \qj-qi\+qj-qb-qi + \qi-qk\-qc+\qh-qk\-$-5e 

c,h,i,j,k,l 

r< r 2qa ~ s - 5e . (3.5) 
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Term 4eD is 



+-q a +qb%acS~hn~ p S~il~ rk 
1 9 9 y(bj)Jin 9 9(kc)Jhl 

£-qa+qb+\qa-qc\+\qh-qn\+\qb-q J \+qj-q t -qn+\qi-qi\+\qk-qc\+qk-qh-qi-$-5< : - 

c,h,i,j,k,lj^h,n^i 

_^ ^2q b -2q d -2q d _ 1 -S-5e ^g 

Term 5D is 

-<; ^-2g m i n {a,c}+2gi-2g min{j , fc }-(5-3e ^g 

The estimates of the remaining terms are obtained as these. The examples include one of 
the terms which limits the estimate for each possible choice of indices a and b. The result 
is, 

s R a b ~^ t 29b ~ 29d ~ 2qd ~ 1 ~ s ~ 5e + ^2 t~ 29inin{a ' c}+29l ~ 29min{j > k ' t ~ s ~ 5e (3 8) 

And 

^2-a a b Sj^a^ _^ ^2g mln{aiij} -2g d -2qr d _i _|_ £-2g max { a ,b}+2gi-2<3 min{j , fe } j^2-o -<5-5e 

^ £2-2q d -2q d _ 1 +2q 1 -(d+7)e ^a-(d+7)e _ Q(+S\ /g g^j 

The estimate of the rest of the terms in / is obtained straightforwardly by checking 
that the exponent of t in each case is strictly positive. The other regularity conditions 
that / should satisfy are shown to hold by equation (31) in [H| and the remarks following 
equation (31). The symmetry of g a b is shown for all efs in subsection |2.5.3j . That the 
Hamiltonian and momentum constraints are satisfied is shown by the direct analogue of 
argument made in section [2^ and the estimate R = o(t~ 2+m ) obtained from equation (|3.8| ). 
The only change is that equation ( |2.80| ) is replaced by 

g ab T c ab C c t- pc r< t~ 2pd - 5 . (3.10) 

To conclude this section we discuss the solution of the velocity-dominated constraints 
for the vacuum equations and D > 4. The case D = 3 could be discussed in a similar way 
but the analogue of the Lichnerowicz equation has a different form and so for brevity that 
case will be omitted. The discussion proceeds in a way which is parallel to that of the last 
section. As already indicated there, the essential task is the analysis of the Lichnerowicz 
equation. In the present case we start with free data g a b, k a b and H where k a b has zero 
divergence. The actual data are defined by g a b = ^ A ^ d ~ 2 ^g a b and k a b = cu~ 2 k a b + Hg a b- 
The constraints will be satisfied is to satisfies the following analogue of the Lichnerowicz 
equation: 

A s co + J0-(-R- g u; + rt.^) - $±Z*Ih*u,& = (3.11) 
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The corresponding equation in the velocity-dominated case is 



(1 — 2 ~ >~ 3d-2 d( d — 2 ) „ 9 d+2 

-k ab k ab u— - — J -H 2 uo^ = (3.12) 



4(d-l) ao 4 

As in the case of Q2.95p it is trivial to solve ( |3.12| ) provided k ab does not vanish at any 
point. To determine the solvability of equation (|3.11|) it is necessary to study the positive 

3d-2 d+2 3d-2 d+2 

zeros of the algebraic expressions x + bx d ~ 2 — ax d - 2 and bx d ~ 2 — ax d ~ 2 where a > and 
b > 0. The second expression is very close to what we had in the velocity-dominated case 
and clearly has a unique positive zero for any values of a and b satisfying the inequalities 
assumed. Looking for positive zeros of the first algebraic expression is equivalent to looking 

d+2 „ 

for positive solutions of x d ~ 2 + ax —6 = 0. Note that the function on the left hand 
side of this equation is evidently decreasing for all positive x, tends to infinity as x — > 
and tends to — b as x — > oo. Hence as long as the constant b is non-zero this function 
has exactly one positive zero, as desired. This is what is needed to obtain an existence 
theorem. It would be desirable to also obtain a uniqueness theorem for the solution of 
( p.ll ). To obtain solutions of the velocity-dominated constraints of the right kind to be 



consistent with Kasner-like behaviour we can use the same approach as in the last section, 
starting with Kasner solutions with an appropriate set of Kasner exponents. 



4 Massless scalar field, D > 3 

Consider Einstein's equations, D > 3, with a massless scalar field as source, the action 



given by SE[g a p] + S<f>[9ap, </>], and d+ 1 decomposition as in section [L6|. The stress-energy 
tensor is 

T, v = ^ V,0 - \g, v ( D V>. (4.1) 

Thus p = \{(d t <P) 2 + o a6 e a (0)e 6 (0)}, j a = -d t 0e a (0), and M\ = g ac e b (<P) e c (0). A 
crucial step in the generalization to arbitrary D > 3 is that the cancellation of terms 
involving d t <f) in the expression for M a b is not particular to D — 4. The scalar field satisfies 
( D ) V a ^ V a 4> = 0, which has d + 1 decomposition 

d 2 t cf> - (trfc)^ = g ab V a V b 0. (4.2) 



Let the Kasner-like evolution equations be equations (|2.7|) - ( |2.9| ), with solution (|2.12|) 
- (|2.14j ) for time coordinate as in section Given a point xq G E, let the neighbourhood 
Uo, the (local) adapted frame and the constant e be as in section |3|. Define °p = \{d t °(f)) 2 
and °j a = —dt°<f)e a ( <j)). The velocity-dominated constraints corresponding to the Hamil- 
tonian and momentum constraints are °C = and °C a = 0, with °C and °C a given by 
equations ( 2.1 8|) and (|2.19| ). For the solution ( p.12 ) - ( 2.14 ) the velocity-dominated Hamil- 



tonian constraint is equivalent to Y^Pa 2 + A 2 = 1. Equations ( |2.21[ ) and ( |2.22| ) are satisfied 
so if the velocity-dominated constraints are satisfied at some t , then they are satisfied 
for all t > 0. The restriction defining the set V is the inequality ( ^.1| ). (If D < 11, then 
satisfying simultaneously (|3TT|), J^Pa = 1 and J^Pa 2 + A 2 = 1 requires A ^ 0. Note that 
conversely, for D = 3, the restrictions defining V are simply equivalent to A ^ 0, since 
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Q3.1D is in this case a consequence of p\ + P2 = 1 and p\+ pi < 1.) The constant a > is 
chosen so that, for all x G U , 4a < 1 + p\ — pd — Pd-i from which it follows that 



a < 2 - 2p d . 



(4.3) 



Now reduce e if necessary so that e < a/ {2d +1). As before, this may in turn require 
shrinking Uq. 

The unknown u = (7%, \ a bc , K a b , ip, u a , x) is related to the Einstein-matter variables 
by ( |2.31|) - ( p.36| ). The quantities A and / appearing in equation ( |1.5| ) are given by the 
evolution equations ( 2.39|) - (|2.43|) and 



t d tX + Px = t a °-\ti k)(A + t?x) + t 2 - p s g ab s V a s V b °<f> + t 2 " c s V a u a . 



(4.4) 



The argument that the matrix A satisfies the appropriate positivity condition is analogous 
to the argument in ||I5 |. Regarding the estimate / = 0(t s ), the estimate t 2- "" 6 s R a b = 
0(t s ) was obtained in equation ( |3.9| ). The estimate t 2 ~ a °' b M a i ) = 0(t 5 ) follows from the 
inequality ( f4.3|) and from qa < Pd- The only other terms in / whose estimates are not 
immediate from the estimates made in [[DJ are the last two terms on the right hand side of 
equation ( |4.4|) . The covariant derivative compatible with the symmetrized metric is used 
in equation 
and (12341 ), 



for convenience. From the estimate s g ab ^ t' 9a 9b ' e p!5| |, equations ( J2.53 ) 

and e c {g [ab) ) ■< t 2q —^- & -\ 



S gab _^ ^-2g min { a , i) }-e 



Therefore, 



S ab S-pc 
9 1 ab 



Sgab S 



9 ch (ea{9(bh)) - le h {g (ab) )) - s g ch f a ah 



-< t 



-2q d -8-3e 



and 



t 2- P SgabS^J^^^ 



J- „h& 



= t 2 -? VHe a (e,(V)) 

^ t 2-2a d -P-S-3e = (t S ), 

= t 2 " C S g ab {e a (w b ) - s T c ab w c } 



ab^c\ 



V)} 



The other regularity conditions that / should satisfy are again shown to hold by equa- 



tion (31) in [15] and the remarks following equation (31). That g ab is symmetric (so that 
equation (|4~4l ) and equation ( [i~2D are equivalent) is shown as in subsection [2.5.3| . That 
the Hamiltonian and momentum constraints are satisfied is shown by the analogue of the 
argument made in section [2^ and the estimate R = o(t~ 2+Vl ) obtained from equation (|3.8| ). 

Note that the case D = 3 of this result has an interesting connection to the Einstein 
vacuum equations in D = 4. As it follows from standard Kaluza-Klein lines, the solutions of 
the latter with polarized U(l) symmetry are equivalent to the Einstein-scalar field system 
in D = 3 (see e.g. [^0| and [41 j, section 5). Hence the result of this section implies that 



we have constructed the most general known class of singular solutions of the Einstein 
vacuum equations in four spacetime dimensions. These spacetimes have one spacelike 
Killing vector. 
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5 Matter fields derived from n-form potentials 
5.1 Equations of motion 

We now turn to the general system , but without the interaction terms "more" . These 
are considered in section [7] below. The action is the sum of (|1.2| ), (|1.3|) and k additional 
terms, each of the form ( |1.4| ). The argument is based on that of section |j. It is enough 
here to note the differences. Furthermore, since there is no coupling between additional 
matter fields, the differences from the argument made in section |3] can be noted for each 
additional matter field independently of the others. Therefore consider the jth additional 
matter field, 

with A an rij-form. This matter field contributes the following additional terms to the 
stress-energy tensor, equation (P~T|), 



f 1 1 

rp _ I J p p ai--a„. p p 



3 ' 



Define £ ai '" an i = yjfg F 0ai '" a ™j e Xj< ^. If rij = 0, £ is a spatial scalar density. Throughout 
this section and the next we use the following conventions. If rij = 0, then P ai ... an . is a 
scalar, g aibl ■ ■ ■ g ajl .b„. = 1, etc. The d + 1 decomposition of the contribution of this matter 
field to the stress-energy tensor is 

„ I J „ „ c«i-«n, ■ cbi—b„ ■ -\j4> 

9 ~ "' + 2g~^. gaibl '" 9an ^ 

+ n , ~ T^ 9 aobo ■ ■■g an ^F ao ... an F bo ... K e x >\ (5.1) 



■>iii , ■ 1 )! ' 



3a = - + — ^F*.. . bn ., (5.2) 
yfgnA 

A/T a — ^ I n i n n n pooi—Cns f>h\—hni 

M b - ( — -.gbhx9c % h^ ■ ■ -gcn-hn^ 3t 3 



g 



3 m 



(d-l) nj \ 

I ( \ „aC r h\l\ r.hn.in. 771 771 

+( — -,9 9 ---g 3 • './', /-,,, bi • 
^ 5\g Coho ■ ■■g" h ' /•;,,..,,. /•'/,,. 7,.,,)' A ''- 



(d-l)(n, + l)! 
The jth matter field satisfies 

(D) V At (F^ 1 " 1 Se A ^) = 0, (5.3) 

(D) V[A...^.] = 0, (5.4) 

with d + 1 decomposition into constraint equations, 

1 n 3 

ea (^-S) + /^^-S+ -£/^ a62 - c -S = 0, (5.5) 

1 i=2 
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e[a(Fb -b nj ]) - - J 2 flabo F \c\bi---b n .] = 0, (5.6) 

and evolution equations, 

d t £ a ^ = -e b (^g bc °g a ^ ■ ■ ■ F Co ... Cn .e x ^) - {Cg bc °g a ^ ■ ■ ■ g a ^ 

+ fbh9 bC0 9 aiCl ■ ■ ■ 9 h(H ■ ■ ■ g^^VdFco...^ e x ^, (5.7) 

d t F ao ... an . = -(n i + l)e [ao (^ aiN -..^ N £ fcl -"S e -^) 



2 ^ -/[ C aoai ^| c || bl |^ 2 |6 2 | ■ - ' A ' rt . (5.8) 

The jth matter field contributes the following terms to the evolution equation (|4.2| ) for </>. 
d 2 t <P - (trk)d t 4> = ... + -^-g aibl ... g b r"'S5^"S e -^ 

XL ^9 aobo ■ ■■g an ^F ao ... an F bo ,e x ^ (5.9) 



2K- + 1)! 

5.2 Velocity-dominated system 

The Kasner-like evolution equations corresponding to this matter field are d t °£ ai ' a " J = 
and d t °F ao ... an , = 0. The quantities °£ ai '" an j and °F ao ... an . are constant in time with 
analytic spatial dependence and both are totally antisymmetric. 

The velocity-dominated matter constraint equations are equations ( |5.5|) and (|5.6| ) with 
°S and °F substituted for £ and F. Since all quantities in the velocity-dominated matter 
constraints are independent of time, if the matter constraints are satisfied at some time 
to > 0, then they are satisfied for all t > 0. This matter field does not contribute to °p. 
Its contribution to °j a is the term shown on the right hand side of equation (|5.2| ) with °g, 
°£ and °F substituted for g, £ and F . The velocity-dominated constraints corresponding 
to the Hamiltonian and momentum constraints are °C = and °C a = 0, with °C and 
°C a given by equations ( |2.18|) and fl2.19|) . Equations (|2.21|) and (|2.22| ) are satisfied, so as 



before, if the velocity-dominated constraints are satisfied at some to > 0, then they are 
satisfied for all t > 0. 

The presence of the matter field puts the following restrictions on the set V JO. 

2p 1 -\ \-2p nj -X j A>0 and 2p 1 -\ h 2p d _ n ._ x + XjA > 0. (5.10) 

The restrictions generalize the inequalities (|2.24j) found for a Maxwell field in 4 dimensions 
and, like them, guarantee that one can asymptotically neglect the p-form A^ in the 
Einstein-dilaton dynamical equations. (For rij = 0, the inequality on the left of ( |5.10| ) is 



— XjA > 0, while for rij = 1 it is 2pi — XjA > 0. For rij — d — 1, the inequality on the right 
is XjA > 0, while for rij — d — 2 it is 2p x + XjA > 0.) 

The constant a is reduced from its value in section |4], if necessary, so that, for all x G Uq, 
o < 2pi + • — h 2p nj — XjA and a < 2p\ + • — h 2pd- nj -i + XjA. If a is reduced, it may be 
necessary to reduce e, and in turn shrink Uq, so that the conditions imposed in section |] 
are still all satisfied. 
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5.3 Fuchsian property - Estimates 

The jth matter field contributes the following components to the unknown u in the Fuch- 
sian equation fll.5Q . 

gai-a nj = . / (5.H) 

Fao—On-A = Fa ---a n . + t (fao—cinj • (5.12) 

Here, j3 = e/100 as above, £ 0l " ,an j is a totally antisymmetric spatial tensor density, so 
contributes Q*) independent components to it, and f> ao - anj is a totally antisymmetric 
spatial tensor, so contributes ( n d + -^j components to u. This is consistent with the form of 
the evolution equations. Note that £ ai "' a "j =0(1) and F ao ... an , = 0(1). 

This matter field contributes additional rows and columns to the matrix A such that 
the only non-vanishing new entries are on the diagonal and strictly positive. Therefore, 
the presence of this matter field does not alter that A satisfies the appropriate positivity 
condition. 

The terms in the source / which must be estimated on account of the jth matter field 
are the following. It contributes terms to the components of / corresponding to K through 
its contribution to MV 

+2-a a b 1 cac2---c„- chi---h n .-Xj(f) 

t -9bhi9c 2 h2 Qcn.hn.^ J t j e 
9 

£— 2 ?max{a,b}+ 2 9a.+ 2 9max{6,fei}H h 2 9majt{c nj . ,h nj y—^jA-ao—nje 

_^ ^2giH \-2q nj -\jA-ao-nje _ Q/^-a -%*e+<j\ _ 0(t & ) (5 13) 

t 2-a« bg ac ghlh . . . gW^ Fchl ... hn F hll ... ln /^ 

_^ ^ 2 ~ 2 9max{a,i)}+ 2 9o— 2 g m in{a,c}— ^minlfe!,^} 2( /min{h„ j } +Aj J 4-«o-(n j +l)e 

_^ £ 2 <?H h2(j d _ nj _i+A JJ 4-ao-(«j+l)e _ Q/^-c«o-(nj+l)e+cr\ _ (J(t & ) (5 14) 

Here it is used that both £ ai " a ™j and F ao ... a ^ are totally antisymmetric, so that the sums 
indicated by a summation symbol are not over all indices. Note that the inequalities ( |5.10| ) 



have been crucially used in getting the estimates (|5.13|) and ( |5.14j) . The desired estimates 
for the other two terms are obtained similarly. 

The terms contributed to the component of / corresponding to \ by the jth matter 
field are obtained by multiplying the right hand side of equation ( p.9[ ) by t 2 ^ 13 . 

t 2 - p -9a lb r--9a n .K.£ n ^ bl '^e- x ^ = 0{t-P- n i* + °) = 0(t s ) (5.15) 

9 3 3 

t 2 -Pg a ^...g a ^F ao ... an F bo ... bn e x ^ = O^"^ 1 ^) = 0(t a ). (5.16) 



3 

The terms in / corresponding to £ ai "' an j for the jth matter field are obtained by multiplying 
the right hand side of equation ( |5.7| ) by t l ~@ . These terms are 0(t~^ 3 ~ & ~^ n3+l ' >e+ ' 7 ) = 0(t s ). 
The terms in / corresponding to (f ao ---a n for the jth matter field are obtained by multiplying 
the right hand side of equation ( |5.8|) by t l ~@ . These terms are 0{t~^~ 5 ~ n: > e+<T ) = 0(t s ). 
Thus the terms which occur in / due to the jth matter field are 0(t 5 ). 
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The time derivative of the matter constraint quantities for the jth field (the left hand 
side of equations ( |5.5|) and ( |5.6|) ) vanishes. If the velocity-dominated matter constraints are 
satisfied, the matter constraint quantities are o(l). A quantity which is both constant in 
time and o(l) must vanish. Therefore the matter constraints for the jth field are satisfied. 

Next the matter terms due to the jth field in the Einstein constraints are estimated, in 
order to verify that they are consistent with equations ( |2.63| ) and (|2.64|) . The contribution 
to the Hamiltonian constraint is, from equation ( |5.1D , 



-g aibl ---9a n . bri £ ai - a ^£ bl -^e- x ^ = 0{r 2 ~ n ^) = o{t- 2+ ^) 1 (5.17) 

g J J 

g aobo ---g an ^F ao ... anj F bo ... bnj e x ^ = 0(r 2 -^ +1 ^) = o(t' 2+ ^). (5.18) 
The contribution to the momentum constraint is 

V9 V 9 J 
+^-(E b ^i F abl ... bn . - °F abv .. bn .) = o(r 1+ ^). 



Estimates related to the determinant which are relevant to ( |5.19| ) are analogues of the 



estimates for d = 3 immediately preceding equation (|2.55|) . The form of these estimates 



for general d will now be presented. These are 1/y/g- l/\f 1Y g = 0(r 1+a °~ dt ), e a (g) 

o(t a °- 5 - de ), 

e a (g) = 0(t 2+a »- 5 - d *), (5.20) 

and 

e a (g- l/2 ) = -01 = 0(r 1+ --^). 
6 Determination of subcritical domain 

The explicit determination of the subcritical range of the dilaton couplings for which the 
inequalities on the Kasner exponents are consistent so that V exists may be a complicated 
matter. We consider in this section a few cases and give some general rules. As in subsection 
2~1, we introduce the metric 



dS 2 = G^dtdtf = £ d Pa 2 - (£ dPaf + (dA) 2 (6.1) 

in the D-dimensional space of the Kasner exponents (p a , A) = (p M ). This metric has again 
Minkowskian signature (— ,+,+,•••,+). The forward light cone is defined by 

twy = o, E^>°- (6-2) 

The Kasner conditions met in the previous section are equivalent to the conditions that 
the Kasner exponents be on the forward light cone (since Pa = 1 can always be achieved 
by positive rescalings). 
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The wall chamber W is now defined by 



Pi < P2 < ■ ■ ■ < Pd (6.3) 
2pi + P2 + ■ • • + Pd-2 > (6.4) 



and, for each p-form, 



pi + p% + ■ ■ • + P nj ~ ^-A > (6.5) 
P1+P2 + ---+ Pd- nj -i + > 0. (6.6) 

These inequalities may not be all independent. The question is to determine the "allowed" 
values of the dilaton couplings for which the wall chamber contains in its interior future- 
directed lightlike vectors. It is clear that this set is non-empty since the inequalities can 
be all fulfilled when the couplings are zero (the p a 's can be chosen to be positive in the 
presence of a dilaton). 

6.1 Einstein-dilaton-Maxwell system in D dimensions 

We consider first the case of a single 1-form in D > 4 dimensions. This case is simple 
because the inequalities ( |6.4|) are then consequences of ( |6.5|) and (|6.6j ), which read 

Pi-y4>0, pi+p 2 + ---+Pd-2 + -A>0. (6.7) 

Furthermore, the number of faces of the wall chamber (defined by these inequalities and 
(|6.3|)) is exactly D and the edge vectors form a basis. Thus, the analysis of subsection |2]4] 
can be repeated. 

A basis of edge vectors can be taken to be 

(0,0, ■•■,0,1,0) (6.8) 

d-k-2 _ d-k-2 2(d-fc-2) k = 12 ... d _ 2{69) 

[ k + l ' ' fc+1 ' ' ' ' X(k+1) h ' ' ' 1 j 

(6.10) 

1 - L ) (6.11) 

In ( |6.9|) , the first k components are equal to — d fc ^ 2 and the next d — k components are 
equal to 1. 

The first vector is lightlike. The k-th vector in the group ( |6.9|) has squared norm 

(d - l)\k 2 - k(d - 3) + d] A(d-k-2) 2 , . 
+ W + l)» • *=W.-,*-* (6-12) 

while (|6.10D and ( |6.11D have norm squared equal to 



-d{d-l) + ^ (6.13) 
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-d(d-l) + 1 " ; , (6.14) 



and 

4(d_ 
~A 2 

respectively. The subcritical values of A must (by definition) be such that at least one 
of the expressions ( |6.12| ), ( |6 . 1 3| ) or fl6.14|) is positive. To determine the boundaries ±A C 



of the subcritical interval, we first note that ( |6.13 ) is positive whenever |A| < Ai, with 



Ai = 2/ \Jd(d — 1). Similarly, (|6.14j ) is positive whenever |A| < A2 with A2 = 2(d — 
2)/yd(d — 1). To analyse the sign of Q6.12 ), we must consider two cases, according to 



whether k 2 — k(d — 3) + d is positive or negative. 

If d < 9, the factor k 2 — k(d — 3) + d is always positive (for any choice of k, k = 
1, 2, • • ■ , d — 2) and the expression (|6.12|) is positive provided |A| < Ilk, with 

= 2(*-*-2) 

'(d - l)[Jfe 2 - ife(d- 3) + d] 



The critical value A c is equal to the largest number among Ai, A2 and H^. This largest 
number is A2 for d = 3,4, 5, 6, TEi for d = 7 and H2 for d — 8. We thus have the following 
list of critical couplings: 

A„ = y|, d = 3 

A c = -|, , = 4 
3 

Ac = 7H' rf = 5 

a c = 1^, d = 6 

x 2v/ 2 w 7 
A c = 7 =-, rf= 7 

4a/2 

A c = ^, d=8. (6.16) 

Note that the value of the dilaton coupling that comes from dimensional reduction of 
vacuum gravity in one dimension higher 



X KK = 2^ (6.17) 

is always strictly greater than the critical value, except for d = 8, where \ KK = A c . (The 
corresponding values of the Kasner exponents are those of the point on the Kasner sphere 
exhibited in |TJ| for D = 10, where all gravitational inequalities are marginally fulfilled.) 



If d > 9, the factor k 2 — k(d — 3) + d is non-positive for 



d - 3 - J(d - 9)(d - 1) d - 3 + J(d - 9)(d - 1) 

^ < fc < ^ (6.18) 
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(this always occurs for k = 3). Thus, the expression ( 6.12| ) is positive for such k's no 
matter what A is. This implies that the critical value of A is infinite, 



A r 



oo, 



d > 9. 



(6.19) 



The fact that D = 10 appears as a critical dimension for the Einstein-dilaton-Maxwell 
system, above which the system is velocity-dominated no matter what the value of the 
dilaton coupling is in the line of the findings of [ 13| , since the edges (|6.9|) differ from those 
of the pure gravity wall chambers only by an additional component along the spacelike 
dilaton direction. 



6.2 Einstein-dilaton system with one ]>form (p ^ 0, p ^ D — 2) 

The same geometrical procedure for determining the critical values of the dilaton couplings 
can be followed when there is only one p-form in the system (p ^ 0, p ^ D — 2), because in 
that case the wall chamber has exactly D faces and the edge vectors form a basis. Indeed, 
the gravitational inequalities (|6.4|) are always consequences of the symmetry inequalities 
( |6.3|) and the form inequalities Q6.5D and (|6.6| ) (for rij ^ and rij ^ D — 2), 

2pi + P2 + ■ ■ • + Pd-2 = (pi + • • • + Pn 3 - y^4) + (pi + Pn,+i + ••• + Pd-2 + yA). (6.20) 

So, if there is only one p-form (with p ^ and p ^ D — 2), the D — 2 symmetry inequalities 
( p.3| ) together with the two form inequalities ( |6.5| ) and ( |6.6| ) completely define the wall 
chamber, which has D faces. We shall not provide an explicit example of a calculation of 
A c for such a system, since it proceeds as for a 1-form. 

When there is more than one exterior form, one can still drop the gravitational inequal- 
ities (if there is at least one p-form with p ^ and p ^ D — 2), but the situation is more 
involved because the inequalities corresponding to different forms are usually independent, 
so that the wall chamber has more than D faces (its intersection with the hyperplane 
J^Pa — 1 is not a simplex). The calculation is then more laborious. The same feature 
arises for a 0-form, which we now examine. 

6.3 0-form in 4 dimensions 

We consider the case of a 0-form in 4 spacetime dimensions. As explained above, we impose 
the condition A ^ to the corresponding dilaton couplingQ Without loss of generality (in 
view of the <fi — ► —<ft symmetry), we can assume A > 0. The inequalities defining the 
subcritical domain relevant to the 0-form case can be brought to the form 

pi > (6.21) 
A > (6.22) 

7 The case A = is clearly in the subcritical region but must be treated separately because there are 
then two dilatons. The Kasner conditions read pi + p2 + ■ ■ ■ + Pd = 1 and p\ + ■ ■ ■ p\ + Af + A| = 1, where 
the scalar fields behave as <p\ ~ A\ hit, cf>2 ~ A2 Int. This allows positive Pi's, which enables one to drop 
spatial derivatives as t — > 0. The system is velocity-dominated. 
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Px+p 2 --A>0 (6.23) 

p 2 - Pi > (6.24) 
P3 ~ P2 > (6.25) 

We denote by a, /3, 7, 6 and e the corresponding border hyperplanes (i.e., a : p\ = 0, 
/3 : A = etc). The inequalities ( |6.21| )-( |6.25| ) guarantee that all potential walls are 
negligible asymptotically. They are independent. The five faces a, (3, 7, S and e intersect 
along the 7 one- dimensional edges generated by the vectors: 

ei = (0,0,l,0)ean/3n7 = an/3n<J = an7n<y = /3n7n<5 (6.26) 

e 2 = (0,1,1,0) G an/3 He (6.27) 
2 

e 3 = (0,1,1, T ) earn 7De (6.28) 
A 

e 4 = (0,0,0, 1) e an (5 He (6.29) 

e 5 = (-1,1,1,0) e /?n 7 ne (6.30) 

e 6 = (1,1,1,0) e p nSne (6.31) 

e 7 = (1,1,1,-) € 7n5ne (6.32) 
A 

Among these vectors, neither e 4 nor e§ bound the subcritical domain since e 4 is such that 
Pi + P2 — (A/2)v4 < (changing its sign would make A < 0), while e$ is such that p± < 
(changing its sign would make P2 — pi < 0). 

The edge-vectors {ei, e 2 , e 3 , e 6 , e 7 } form a complete (but not linearly independent) set. 
Any vector can be expanded as 

v = v x ei + v 2 e 2 + v 3 e 3 + v 6 e 6 + v 7 e 7 (6.33) 

The coefficients v±, v 2 , v 3 , v e , Vj are not independent but can be changed as 

v 2 — > v 2 + 2k, v 3 v 3 — 2k, Vq v 6 — k, v 7 —> v 7 + k (6.34) 

For (|6.33| ) to be interior to the wall chamber, the coefficients V\, v 2 , v 3 , v 6 and v 7 must 
fulfill 

vx > 0, v 2 + v 3 > 0, v 2 + 2v 6 > 0, v 3 + 2v 7 > 0, v e + v 7 > 0. (6.35) 

Using the above redefinitions, which leave the inequalities invariant, we can make v a > 0, 
A = 1,2,3,6,7, with at most two tu's equal to zero. Indeed, let s = min(f 2 , v 3 , 2f 6 , 2v 7 ). 
Assume for definiteness that s = v 2 (the other cases are treated in exactly the same way). 
One has then t> 2 < 2v 7 . Take 2k = —s in the redefinitions ( |6.34|) . This makes f 2 equal 



to zero and makes v 7 equal to v 7 — (v 2 /2) > 0. Because of (|6.35|) , the new v 3 and v Q are 



strictly positive, as claimed. Thus, one sees that any vector in the wall chamber can be 



expanded as in (|6.33| ) with non- negative coefficients. But the vectors e\, e 2 , e 3 , e 6 and e 7 



are all future-pointing and timelike or null when A > y8/3. It follows that for such A's, 
there is no lightlike direction in the interior of the wall chamber. Conversely, if A < ^8/3, 
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the vector e-j is spacelike and one can find an interior vector ae\ + [3e2 + t-j (a, (3 > 0) that 
is lightlike. We can thus conclude: 



for a 0-form in 4 dimensions, (6.36) 



i.e., the system is velocity-dominated for |A| < y8/3. 

The action for the matter fields in the case of a 0-form A coupled to a dilaton </> is 

S4g al 3, 0, A \ = -\j °y + e "% A 9^)^=9 d 4 x (6.37) 

Note that this is the action for a wave map (also known as a nonlinear cr-model or hyper- 
bolic harmonic map) with values in a two-dimensional Riemannian manifold of constant 
negative curvature. Its curvature is proportional to A 2 . Thus we obtain an interesting 
statement on velocity-dominated behaviour for the Einstein equations coupled to certain 
wave maps. Note for comparison that wave maps in flat space occurring naturally in 
the context of solutions of the vacuum Einstein equations with symmetry, for instance in 
Gowdy spacetimes (cf. p4]|), are defined by a Lagrangian of the above type (using the flat 
metric) with A = 2. 

6.4 Collection of 1-forms 

We now turn to a system of several 1-forms. It is clear that if these have all the same 
dilaton coupling, as in the Yang-Mills action (|1.7| ), then, the critical value of A is just that 
computed in ( |6.16|) and ( |6.19|) since each form brings in the same walls. The situation is 



more complicated if the dilaton couplings are different. One could naively think that the 
subcritical domain is then just the Cartesian product of the individual subcritical intervals 
[— \v\\v>], but this is not true because the intersection of the wall chambers associated 
with each 1-form may have no interior lightlike direction, even if each wall chamber has 
some. 

This is best seen on the example of two 1-forms in D spacetime dimensions with opposite 
dilaton couplings. The relevant inequalities, from which all others follow, are in this case 

Pi - > 0, Pl + ^A>0 (6.38) 

Pi<p 2 < ■■■ <Pd (6.39) 

and can be easily analysed because they determine, in this particular instance, a simplex in 
the hyperplane Y,Pa — 1- ft follows from ( p. 38] ) that p\ > 0. The edge- vectors can be taken 
to be (0, • ■ ■ , 0, 1, ■ • • , 1, 0) {k zeros, d — k ones, k — 1, • • • , d — 1) and (1, 1, • • ■ , 1, ±2/A). 
The first d — 1 edge- vectors are timelike or null, while the last two are spacelike provided 
-d(d - 1)A 2 + 4 > 0. This yields 

2 

A c = = for two 1-forms with opposite dilaton couplings (6.40) 

y/d{d-l) 

Accordingly, A c is finite for any spacetime dimension (and in fact, tends to zero as d — > oo), 
even though A c = oo for a single 1-form whenever d > 8. 
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7 Coupling between the matter fields 



The actions for the bosonic sectors of the low-energy limits of superstring theories or M- 
theory contain coupling terms between the p-forms, indicated by "more" in These 
coupling terms are of the Chern-Simons or the Chapline-Manton type. In this section, we 
show that these terms are consistent with the results obtained in section |j| in that they 
are also asymptotically negligible in the dynamical equations of motion when the Kasner 
exponents are subject to the above inequalities (|6.3|) - (|6.6|) . 

More precisely, the form of the velocity-dominated evolution equations and solutions 
are in each case exactly as in section [5[ The velocity-dominated matter constraints have 
additional terms, but as before, the velocity-dominated matter variables (besides the dila- 
ton) are constant in time, so if the constraints are satisfied at some t > they are satisfied 
for all t > 0. The quantities °p and °j a are defined exactly as in section |^. Since the 
velocity-dominated evolution equations are also the same, there is nothing additional to 
check concerning the velocity-dominated Hamiltonian and momentum constraints. 

Turning now to the exact equations, the restrictions defining the set V are unchanged 
from section |5[ The form of the evolution equation for the dilaton is unchanged. The form 
of the stress-energy tensor is also unchanged, and so the form of the Einstein evolution 
equations and the Einstein constraints is unchanged. The additional matter field variables 
considered in section [5| are still all 0(1), so estimates of terms involving the matter fields 
do not need to be reconsidered, as long as their form has not changed, for instance, in 
the argument that the Einstein constraints are satisfied. That the matter constraints are 
satisfied follows as in the other cases, once it is verified that their time derivative vanishes 
and that they are o(l). Since so much of the argument is identical to that of section |||, we 
only point out the few places where there are differences. 

7.1 Chern-Simons terms 

First we consider the coupling of i of the additional matter fields via a Chern-Simons term 
in the action. These additional matter fields should be such that 

i-1 + = (7.1) 

3=1 

The Chern-Simons term which is added to the action is 

5 cs [4?.. 7ni , • ■ ■ , 4g. = J A« A dA^ A • ■ ■ A <L4». (7.2) 

The variation of this term with respect to both the metric and the dilaton field, 0, vanishes. 
The matter equation ( |5.4| ) is unchanged, since it is still the case that = dA^ for all 
j. But the equation ( |5.3|) for each of the i coupled matter fields acquires a non- vanishing 
right hand side. 

(jD) V M (F (i)M ' /1 '"^e A ^) v /3 ^ = ()<"'' " /-' :|: ■ • . FV~V FV +1 ') ■ ■ -F® (7.3) 
Here e° = 1 and Cj is a numerical factor. Next, considering the d + 1 decomposition of 



equation ( |7.3|) , the constraint equation (|5.5[) , for the jth coupled matter field, acquires the 
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following term on its right hand side, 



- r> • • • F ij ~ 1} F ij+1) ■ ■ ■ F {i) (7.4) 

Here all indices which are not explicit are spatial. So, only magnetic fields appear in (|7.4j ). 
The following term is added to the right hand side of the evolution equation ( |5.7| ) for the 
jth coupled matter field. 

-C 3 { (n m + i) £ -*'-^-«i-S-F (1) ■ • ■ 

m=l 

([,{ ("»> ■ 1) ( 11 ■■F;; ) ...F; J (7.5) 

m=j+l 

Again, all indices which are not explicit are spatial. There is in each term only one electric 
field. 

The velocity-dominated matter constraint equations for the jth coupled matter field can 
be obtained from the "full" matter constraint equations for the same field by substituting 
the velocity-dominated quantities for all variables. 

The only additional terms occurring in / are due to equation ( [7.5| ). The form of the 
mth term on the right hand side of equation ( |7.5| ) is just like the form of the terms on 
the right hand side of equation ( |5.8| ) for the mth coupled field. The factors which differ, 
comparing the mth term of ( [7.5|) to equation ( |5.8| ) for the mth field, are 0(1). Since in both 
cases a factor of t 1 13 is added in order to obtain the terms appearing in /, the estimate 
that the additional terms in / due to the Chern-Simons coupling are 0(t s ) is obtained just 
as the corresponding previously obtained estimates. 

7.2 Chapline-Manton couplings 

Next we consider Chapline-Manton couplings. For definiteness, we treat two explicit ex- 
amples, leaving to the reader the task of checking that the general case works in exactly 
the same way. The first coupling is between an n-form A and an (n + l)-form B and is 
equivalent to making B massive. Let F = dA+B and H = dB. The gauge transformations 
are B B + dr), for arbitrary n-form rj, and A — > A — rj + d'y, for arbitrary (n — l)-form 7. 
(If n = 0, then d'y is replaced by a constant scalar and we require that the corresponding 
constant, A^, in the coupling to the dilaton be nonzero.) The form of the action is the 
same as in section |5|, but since F now depends on B and not just on A, the variation of 
the action with respect to B acquires an additional term. Also, it is now the case that 
dF = H. 

The matter equation ( |5.3|) is unchanged for F and equation ( |5.4|) is unchanged for H. 
Equation ( |5.3| ) for H and equation ( |5.4j ) for F are now as follows. 
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Define £ a i- a ™ = F * 1 '"^*** and £> a °- a ™ = ^/g H 0a °- •°»e As ^. The matter constraint 
equations which are affected are 

1 n 

e[a(i<bo-b„]) — 2^ — ~/[a6 -^c|bi-6 n ] = ~ 2#ab -b„, (7-9) 



The additional term which appears on the right hand side of equation (57?) for "D 

Jgcf° h °...cf"->»F ho ... hn e x *+. (7.10) 



The additional term which appears on the right hand side of equation ( |5.8| ) for F ao ... an is 

^g aobo ---g anb V b °- b "e-^. (7.11) 

The velocity-dominated matter constraint equations which are affected can be obtained 
from equation (|7.8| ) and ( [7.9|) by substituting the corresponding velocity-dominated quan- 
tities for all variables. The only additional terms occurring in / are due to equations ( |7.10| ) 



for T> and ( [7.11 ) for F. The form of the additional terms in these equations is just like 



the form of the terms which appear in equations ( |5.7| ) for 8 and in ( |5.8| ) for H. Therefore 
the estimate that the additional terms are 0{t & ) is obtained just as the corresponding 
previously obtained estimates. 

The second Chapline-Manton type coupling is between an n-form A and a (2n)-form 
B. Let F = dA and H = dB + A A F. The gauge transformations are A — > A + d'y, for 
arbitrary (n — l)-form 7, and B — ► B + dr) — 7 A F, for arbitrary (2n — l)-form rj. (If n = 
the gauge transformations are A —> A + C and B — > B + D — CA for constant scalars C and 
.D and we require both 7^ and also Xb 7^ 0.) The form of the action is again the same 
as in section §. Define £ a ^ a ™ = Jg F^-^e^ and V a ^ a ^ = ^/g H 0ai - a ^ e XB ^. The 
matter equations ( [5.3D for F and Q5.4| ) for i7 are affected, only if n is odd. The equation 
for F which is affected (if n is odd) and its d + 1 decomposition are 

9 

V (F t±Ul '" Un e^ A ^) = -gliv\—VnO\—an jp e Xs ^ (7 12) 

(n + 1)! 

1 " 2 
ea(r b 2 ...b ?l) + f C aS ab 2 - bn + J2f b a £ ab2 - C " bn = ——V ab ^ h ^ h "F ahl ... hnl (7.13) 

l i=2 [n + I). 

d t S n = ^ pai-a»b r -.b n £ d -c n (7M) 

^ ^9 bh0 g aihl ■ ■ -g^g^+i ■ ■■g Cnh2n H h0 ... h2n F bci ... Cn e XB t 



(n + l)! 
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The equation for H which is affected (if n is odd) and its d + 1 decomposition are 



(D) V tt (2n + l)! 

(2n + 1) (2n + l)! 

e[a(#b -&2„]J ^ /[afeo^klbi-fen] - ( n _|_ _|_ ]\[^[»&o-6„-l^6 n -62„]) I'-J-DJ 

Q rr _i (2n + 2)! cb\---b n —XaQt? (7 M\ 

u t rL a ---a 2n -Jgijl + + 1)! t a °' bl ' ' ' '^ an - 1 \ bn \ r a n -a 2 n}- \ ' ■ 1 ' ) 

The velocity-dominated matter constraint equations which are affected can be obtained 
from equation (|7.13|) and ( [7.1 6|) by substituting the corresponding velocity-dominated 
quantities for all variables. The only additional terms occurring in / are due to equa- 
tions (|7. 14j) for £ and (|7.17|) for H. Here again, the estimate that the additional terms 
in / are 0(t s ), is just as the estimate of terms appearing already in section BL either in 
equation (|5.7|) for V or in equation (|5.8|) for F. 



8 Yang-Mills 

We complete our analysis by proving that Yang-Mills couplings also enjoy the property 
of not modifying the conclusions. The action is ( |1 . 7| ) , with a Yang-Mills field as source 
in addition to the scalar field considered in section [| and with |A| < A c . The argument 
is again based on that of sections Q - |5] and it is enough here to note differences. The 
main one is that one must work with the vector potential instead of the fields themselves, 
because bare A's appear in the equations. We could, in fact, have developed the entire 
previous analysis in terms of the vector potentials, thereby reducing the number of matter 
constraint equations. We followed a manifestly gauge-invariant approach for easing the 
physical understanding, but this was not mandatory. The stress-energy tensor is 

= {D) V M ^ V v <f> - \g,u {D) V a (D) V> + [F m ■ F v a - ^F af3 ■ F^]e x t (8.1) 
We work in the temporal gauge, A = 0. The matter fields satisfy the following equations. 

(D) V Q (D) V> - ^F af3 ■ F a ^e x<t> = (8.2) 

(D) V M (F^e^) + [Ap, F^\e M = 0, (8.3) 

F^ u = d^Av - d u A /I + [An, A,]. (8.4) 

The Lie Bracket has no intrinsic time dependence. The d+1 decomposition of the stress- 
energy tensor is expressed in terms of the spatial tensor density S a = y/g F 0a e x< ^ and the 
antisymmetric spatial tensor F^. 

P = ^{(ft^) 2 + ^e a (^)e^) + - 9at r • ^e- A ^ + ^V"^, • F 6A e^}, (8.5) 
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3a 



-d t <f)eM + ^-£ b -F ai 



g ac e b (<P) e c (0) 



-{g bc £ a • £ c - b\g ch 8 c ■ £ h }e~ x ^ 
9 A 



+{g ac g M Fch ■ F bl - -5 a b g ch g^F a ■ F hj }e x t 

The matter constraint equation is 

e a (£ a ) + f b ba £ a +[A a ,£ a } = 0. 
The matter evolution equations are 



(8.6) 



(8.7) 



d 2 (f> - (trk)d t 



2° 6 V O V 6 + ^g ab £ a ■ £ b e- x * - \g ab g ch F ac ■ F bh e 
2g 4 



\(j> 



d t £ a 



d t A a 



e b (y/gg°*g bh F ch e x *) + {flg ac + ^g ic )Vgg bh F ch e^ 
-^g ab £ b e- x t 

\fg 



(8.9) 
(8.10) 
(8.11) 



Note that we use as basic matter variables A a and £ b (the quantity F ab being then defined 

in terms of A a as F ab = d a A b - d b A a + [A a , A b }). 

The Kasner-like evolution equations are equations (|2.7|) - (|2.10|) and d t °A a = 0. We 
consider analytic solutions of the Kasner-like evolution equations of the form ( 2.12| ) - (|2.15|) 
along with the quantity °A a which is constant in time. Given a point xq G S, we use an 
adapted spatial frame on a neighbourhood of xo, Uq, as in section |3|. Thus, °g a b(to) and 
K a b are specialized as in that section. There is one velocity-dominated matter constraint 
equation, obtained from equation (|8.8[) by replacing £ a and A a with °£ a and °A a . If the 
velocity-dominated matter constraint is satisfied at some time t > 0, then it is satisfied 
for all t > 0. Define 



V = 7;(dt°<f>Y 



J a 



e a ( 



+ 



Ocb Op 



(8.12) 
(8.13) 



The velocity-dominated Einstein constraints are defined as in the other cases. Equa- 
tions ( |2.21| ) and ( |2.22| ) are again satisfied, so if the velocity-dominated constraints are 
satisfied at some to, then they are satisfied for all t > 0. The restrictions defining the set 
V are as in section |5], with rij = 1. The relation of the unknown, u, in equation ( |1.5| ) to 
the Einstein-matter variables is given by equations (|2.31| ) - ( [2. 371) and 

A a = °A a + tP<p a . (8.14) 

The quantities A and / in equation (|1.5| ) are given by equations (|2.39|) - ( [2.43j ) and 

td t x + Px = t a °- /3 (tTK)(A + t l3 x)+t 2 - f}S g abS V a s V b °^ + t 2 ^ s V a u a 
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+t2 ~ P[ Yg 9abSa ' ~ A 9abgChFac ' Fbhe ^ } ' (8 ' 15) 

+Ul9 ac + \f&ftV9ST*Fche»}, (8.16) 
td t <p a + P<Pa = -t l -?^-g ab £ b e- x t (8.17) 

The estimate that / = 0(t s ) is obtained as before, using £ a = 0(1) and F ab = 0(1). 
The matter constraint quantity, the left hand side of equation ( |S.8|) , is o(l) and its time 
derivative vanishes, so the matter constraint is satisfied. The estimate of the matter terms 
in the Einstein constraints is obtained as in section |5| for rij — 1. 

To conclude: the whole analysis goes through even in the presence of the Yang-Mills 
coupling terms and the system is asymptotically Kasner-like provided |A| < A c , where A c 
is the same as in the abelian case and explicitly given by fl6.16|) and (|6.19f) . 



9 Self-interacting scalar field 

Consider Einstein's equations, D > 3, with sources as in sections |], |5|, [7| or [| except 
that the massless scalar field, 0, is replaced by a self-interacting scalar field. That is, 
the expression (|1.8| ) is added to the action. Solutions with a monotone singularity can 
be constructed as in sections [| - BL with assumptions regarding the function V{4>) which 
appears in ( |1.8|) given below. There is no change in the velocity-dominated evolution 
equations and solutions, nor in the velocity-dominated constraints. The only change to 
equation (|1.5|) is that two new terms appear in /. There is a new term, t 2 ~ a °6 a b 2 V(cj>) /{D— 
2), on the right hand side of the evolution equation for K, a t, (through M a f,). There is also 
a new term, — t 2-/3 V(0), on the right hand side of the evolution equation for x- F° r 



equation ([L5|) to be Fuchsian, it must be the case that / = 0(t ) and, in addition, that / 



satisfy other regularity conditions [ITR EH] . 



Some examples were considered in A trivial example is obtained by taking V 



to be a constant. Then the equation for the scalar field is not changed by the potential 
while its effect on the Einstein equations is equivalent to the addition of a cosmological 
constant. Thus we see that the analysis of [|15j generalizes directly to the case of the 
Einstein-scalar field system with non-zero cosmological constant. Of course the analogous 
statement applies to the other dimensions and matter fields considered in previous sections. 
To get another simple example take V(<f>) = \<p p for a constant A and an integer p > 2. 
Showing that the equation is Fuchsian involves examining the expression 

V(Alnt + B + t p i/j) = X(A\nt + B + t^ipy (9.1) 

and corresponding expressions for the first and second derivatives of V. Of course in this 
particular case these are given by multiples of smaller powers of t. The aim is to estimate 
these quantities by suitable powers of t. In this case a Fuchsian system is always obtained. 
A linear massive scalar field is obtained by choosing p = 2. Another interesting possibility 
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is to choose V(4>) = e A< ^ for a constant A, in which case the derivatives of V are also 
exponentials. Then 

V(A\nt + B + t^) = e XB t XA exp(Xt^) (9.2) 

Note that such an exponential potential can be (formally) generated by adding, as matter 
field, a d— form A^...^ with dilaton coupling A^ = —A. Indeed, eliminating the field- 
strength F = dA (which satisfies e Xd ^F = Cr], where C is a constant and r\ the volume 
form), leads to a term in the action proportional to e~ Xd< ^C 2 . A Fuchsian system is obtained 
provided the general "electric" p— form condition ( 5.10Q (with rij = d), 2p% + • • • + 2pd — 



XdA > is satisfied, i.e. (after using pi + ■ ■ • + Pd — 1 and A^ = —A) provided XA > —2. 
This therefore yields a restriction on the data. 

More generally, it is enough to have a function V on the real line which has an analytic 
continuation to the whole complex plane and which satisfies estimates of the form 

t 2 - Cl V(A\nt + B + t^) = 0(1), 

t 2 ~ C2 V\A\xit + B + t^) = 0(1), (9.3) 
t 2 - C3 V"(Alnt + B + t /3 i>) = 0(1), 

for some positive numbers c\, C2 and C3. Here A and B are the analytic continuations of 
A(x) and B(x), to some (small, simply connected) complex neighbourhood of the range of 
a coordinate chart. And ip lies in some region of the complex plane containing the origin. 
For / to be regular, it must be the case that ci > oq and C2 > f3, which can be achieved by 
reducing e, if necessary, and also possibly Uq, so that previous assumptions are satisfied. 
By taking suitable account of the domains of the functions involved it is also possible to 
obtain an analogue of this result when the functions V and V are only defined on some 
open subsets of R and C. 

The only other change to the construction given in sections f|- § is that p — > p + V(4>). 
It is still the case that f iT fJ,v = 0, so equations (|2.59|) and (|2.60|) are satisfied. Equation 
( p.63| ) is satisfied due to the assumptions concerning V(4>), so the Einstein constraints are 
satisfied. 



10 Conclusions 

Our paper establishes the Kasner-like behaviour for vacuum gravity in spacetime dimen- 
sions greater than or equal to 11, as well as the Kasner-like behaviour for the Einstein- 
dilaton-matter systems with subcritical dilaton couplings. Our results can be summarized 
as follows 

Theorem 10.1 Let S be a d-dimensional analytic manifold, d > 10 and let (°g a b, °k a i > ) 
be a C w solution of the Kasner-like vacuum Einstein equations on (0, 00) x S such that 
ttr°k = —1 and such that the ordered eigenvalues of —t°k ab satisfy l+pi — Pd~ Pd-i > 0. 

Then there exists an open neighbourhood U of {0} x E in [0, 00) x S and a C w solution 
(g a b,k a b) of the Einstein vacuum field equations on U H ((0, 00) x E) such that for each 
compact subset K C S there are positive real numbers a a b for which the following estimates 
hold uniformly on K: 
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1. °g ac g cb = n + 0(^0 

2. k\ = °k a b + o{t- 1+aa b) 

Theorem 10.2 Let E be a d- dimensional analytic manifold, d > 2 and let 

X = (°g a b, °k ab , °0, gv) ai '" an 3 ( ); with j taking on values 1 through k for some 

nonnegative integer k (possibly 0, in which case j takes on no values), < rij < d— 1. Let 
Xj be constants in the subcritical range. Let X be a C w solution of the Kasner-like Einstein- 
matter equations on (0, 00) x X such that t tr°k = — 1 and such that the ordered eigenvalues 

of -t °k ab satisfy l+p 1 -p d - p d ^ > and, for each j, 2pi H h 2p n . - Xj t d t °<f) > 

and 2px + '~ + 2Pd- nj -i + Xjtd^cf) > 0. 

Then there exists an open neighbourhood U of {0} x E in [0, 00) x S and a C w solution 
{g a b-, k a b, 0, £ ( -^ ai "' a "j ', F^}..^ ) of the Einstein-matter field equations on U H ((0, 00) x S) 
suc/i i/iai /or eac/i compact subset K C S i/iere are positive real numbers (3, a a b , with 
(3 < a a b , for which the following estimates hold uniformly on K: 

1. °g ac g C b = 6% + o(t aa b ) 

2. k\ = °k a b + o(t- 1+aa b) 

3. = °0 + o{tP) 

Remarks 

1. Corresponding estimates hold for certain first order derivatives of the basic unknowns 
in Theorems |10.1| and |10.2| (cf. Theorem 2.1 in [IBj]). These are the derivatives which 



arise in the definition of new unknowns when second order equations are reduced to 
first order so as to produce a first order Fuchsian system. 

2. Our analysis shows that a solution of the full subcritical Einstein-matter equations 
satisfying the estimates given in the theorems and the corresponding estimates for 
first order derivatives just mentioned is uniquely determined by the solution of the 
velocity-dominated equations (the integration functions are included in the zeroth 
order, Kasner-like solutions; the deviation from them is uniquely determined). 

3. The Einstein-matter field equations may include interaction terms of Chern-Simons, 
Chapline-Manton and Yang-Mills type, and the scalar field may be self-interacting, 
with assumptions on V(<f>) as stated in section |^. If the jth field is a Yang-Mills 
field, then F^ b is obtained from A$ and °F a J b is obtained from through equa- 
tion ( |8.4| ). Note that the condition on tr°A; which is assumed in both theorems can 
always be arranged by means of a time translation. 
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4. The spacetimes of the class whose existence is established by these theorems have 
the desirable property that it is possible to determine the detailed nature of their 
singularities by algebraic calculations. This allows them to be checked for consistency 
with the cosmic censorship hypothesis. What should be done from this point of view 
is to check that some invariantly defined physical quantity is unbounded as the 
singularity at t — is approached. This shows that t = is a genuine spacetime 
singularity beyond which no regular extension of the spacetime is possible. For this 
purpose it is common to examine curvature invariants but in fact it is just as good 
if an invariant of the matter fields can be found which is unbounded in the approach 
to t — 0. This is particularly convenient in the cases where a dilaton is present. 
Then V a 0V a is equal in leading order to the corresponding velocity-dominated 
quantity and the latter is easily seen to diverge like t -4 for t — > 0. The vacuum case 
is more difficult. It will be shown below that the approximation of the full solution 
by the velocity-dominated solution is sufficiently good that it is enough to do the 
calculation for the velocity-dominated metric. This means that it is enough to do 
the calculation for the Kasner metric in D dimensions. Note that the Kasner metric 
is invariant under reflection in each of the spatial coordinates. Hence curvature 
components of the form Ro abc vanish, as do components of the form i?oa06 with a ^ b. 
Hence the Kretschmann scalar R al3 " f6 Ra^s is a sum of non-negative terms of the form 
R abcd R abcd and (i? a 0a o) 2 - In order to show that the Kretschmann scalar is unbounded 
it is enough to show that one of these terms is unbounded. A simple calculation 
shows that (-R a o a o) 2 = P 2 (l _ Va) 2 t~ A in a Kasner spacetime. Thus the curvature 
invariant under consideration can only be bounded as t — > if all Kasner exponents 
are zero or one, which does not occur for the solutions we construct. To see that 
the approximation of the full solution by the velocity-dominated solution is valid for 
determining the asymptotics of the Kretschmann scalar it is enough to note that all 
terms appearing in the Kretschmann scalar which were not just considered are o(t~ 4 ). 
Only two estimates additional to those already obtained are needed - for these, the 
estimates R abc h = 0{t~ 2+t ) and V a k b = 0(t~ 2+e ) are sufficient. Both of these 
estimates are straightforward to obtain. The main input is Y c ah = o(t~ 1+4a ~ 2e ~ s ) 
(i.e., the connection coefficients do not need to be expanded). The expression for the 
Kretschmann scalar is 

4((tr k)k\ - k a c K c b ){{trk)k b a - k b h k\) + (k\k\ - k a h k c b )(k b a k h c - k\k b c ) 
+4{(R\ - M\)(R b a - M b a ) + 2{R\ - M a b )((trk)k b a - k b h k\) 
-2{V a k b c ){V h k\)~g ah ~ 2{V a k b c ){V b k\)~g ch } - ~g ab ~9 ch R aa i R bhj i 
+2R abc \k%k b h -k\k b t )g cl - 

Apart from the Kasner terms (which can each be written as two factors, with each 
factor 0(t~ 2 )), the remaining terms can each be written as two factors, with each 
factor 0(t~ 2 ) and at least one of the two factors o(t~~ 2 ). 

5. We have constructed large classes of solutions of the Einstein-matter equations with 
velocity-dominated singularities for matter models defined by those field theories 
where the BKL picture predicts that solutions of this kind should exist. No sym- 
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metry assumptions were made. When symmetry assumptions are made there are 
more possibilities of finding specialized classes of spacetimes with velocity-dominated 
singularities. See for instance where there are results for the Einstein-Maxwell- 
dilaton and other systems under symmetry assumptions. There are also results for 
the case where the Einstein equations are coupled to phenomenological matter mod- 
els such as a perfect fluid and certain symmetry assumptions are made. For one of 
the most general results of this kind so far see . 



When solutions are constructed by Fuchsian methods as is done is this paper there 
is the possibility of algorithmically constructing an expansion of the solution about 
the singularity to all orders which is convergent when the input data are analytic, 
as in this paper. (If the input data are only C°° the expansion is asymptotic in a 
rigorous sense when Fuchsian techniques can be applied.) At the same time, there 
is the possibility of providing a rigorous confirmation of the reliability of existing 
expansions such as those of p4j and [25]. This is worked out for the case of [24| in 

Pi- 
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